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^"^ ' Abstract. In this paper wc study the convergence of an iterative algorithm 

^ ^ ' for finding zeros with constraints for not necessarily monotone set- valued oper- 

^ ^ ators in a reflexive Banach space. This algorithm, which we call the proximal- 

Q , projection method is, essentially, a flxed point procedure and our convergence 

•^^ ' results are based on new generalizations of Lemma Opial. We show how the 

~ ' ' proximal-projection method can be applied for solving ill-posed variational in- 

1/^ ^ equalities and convex optimization problems with data given or computable by 

approximations only. The convergence properties of the proximal-projection 
method we establish also allow us to prove that the proximal point method 
(with Bregman distances), whose convergence was known to happen for max- 
T-f\ , imal monotone operators, still converges when the operator involved in it is 

, , monotone with sequentially weakly closed graph. 
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1. Introduction 



In that follows X denotes a real reflexive Banach space with norm || • || and X* de- 
notes the (topological) dual oiX with the dual norm ||-||^. Let f : X -^ (—00, -l-oo] 
V/^ , be a proper, lower semicontinuous convex function with domain dom/. Then, the 

CO ' Fenchel conjugate /* : X* — > (—00, -l-oo] is also a proper lower semicontinuous 

^^ , convex function and /** :— (/*)* = /. We assume that / is a Legendre function in 

the sense given to this term in [19l Definition 5.2], that is, / is essentially smooth 
and essentially strictly convex. Then, according to |19l Theorem 5.4], the func- 
tion /* is a Legendre function too. Moreover, by [T^ Theorem 5.6 and Theorem 
5.10], both functions / and /* have domains with nonempty interior, are (Gateaux) 
differentiable on the interiors of their respective domains. 



(1.1) ranV/ = domV/* = intdom/* = dom 9/* 
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(1.2) ranV/* = domV/ = intdom/ = domdf, 

and 

(1.3) v/ = (vr)-i. 

With the function / we associate the function Wj : X* xX ^ (— oo, +00] defined 

by 

(1-4) Wf{^,x) = f{x)-{t^x}+r{0- 

By the Young- Fcnchel inequahty the function Wf is nonnegative and domWf = 
(dom/*) X (dom/). It is known (see [1], [2], [18], [36j and see also Section 2 
below) that, for any nonempty closed convex set E contained in X such that E n 
intdom/ 7^ 0, the function Proj^ : intdom/* -^ X given by 

(1.5) Proj;^^ = argmin{W^/(^,a;) : x e E} , 

is well defined and its range is contained in .E n intdom/. In fact, this function is 
a particular proximal projection in the sense given to this term in [20j which was 
termed in [36] projection onto E relative to / because in the particular case when 
X is a Hilbert space and f{x) = ^ ||x|| the vector Proj;^^ coincides with the usual 
(metric) projection of ^ onto E. 

In this paper we are interested in the following problem: 

Problem 1.1: Given an operator A : X ^ 2^ and a nonempty closed subset 
C of X such that 

(1.6) 7^ CndomA C intdom/, 

find X € C such that 0* G Ax, where 0* denotes the null vector in X*. 

Our purpose is to discover sufficient conditions for the following iterative procedure, 
which we call the proximal-projection method, 

(1.7) x° e Co n dom A n int dom / and 

x'^+i e Proj^,+,ndo„.A(V/(^'^)-^:^'), VfceN, 

to generate approximations of solutions to the Problem 1.1 when {Cfc}j.gpj is a 
sequence of closed convex subsets of X approximating weakly (see Definition 4.2 
below) the set C under the following conditions of compatibility of / and Ck with 
the data of Problem 1.1: 

Assumption 1.1: For each k e N, the set Ck H dom A is convex and closed, 

(1.8) CCCfc and (V/ - A)(Cfe) C intdom/*. 



In this paper (II. 6|) and Assumption 1.1 are standing assumptions, even if not explic- 
itly mentioned, whenever we refer to the Problem 1.1 or to the proximal- projection 
method. In view of (|1.6p and (|1.8p the sets Ck H dom^ n intdom/ are necessarily 
nonempty and, consequently, (V/ — A)(Cfe) is nonempty too. This fact. Assumption 
1.1 and Lemma 2.1 below which ensures that 

domProj^^^j^j^^^ = intdom/* and ranProj^^^^^,^^^ C Cfc H dom^ n intdom/, 

taken together, guarantee that the procedure of generating sequences in (I1.7P is 
well defined. 

The proximal-projection method described above is a natural generalization of 
the method of finding zeros of linear operators due to Landweber [IS] , of Shor's [7D] 
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(see also [71]) and Ermoliev's [44] "gradient descent" methods for finding uncon- 
strained minima of convex functions and of the "projected-subgradient" method 
for finding constrained minima of convex functions studied by Polyak '62] . These 
methods inspired the construction of a plethora of algorithms for finding zeros of 
various operators as well as for other purposes. Among them are the algorithms 
presented in [1], [3], [4], i, [7], [S], [9], [lO], [H], [E], [13], [H], [TT, [20;, [29], 
[35j . [36], [69j which, in turn, inspired this research. The main formal differences 
between the proximal-projection method (jl.7p and its already classical counterparts 
developed in the 50-ies and 60-ies consist of the use of the proximal projections in- 
stead of metric projections and of projecting not on the set C involved in Problem 
1.1, but on some convex approximations Ck of it. The use of proximal projec- 
tions instead of metric projections is mostly due to the fact that metric projections 
in Banach spaces which are not Hilbcrtian do not have many of those properties 
(like single valuedness and nonexpansivity) which make them so useful in a Hilbert 
space setting for establishing convergence of algorithms based on them. As far as 
we know, the idea of using proximal projections instead of metric projections in 
projected-subgradient type algorithms goes back to Alber's works [1], [2]. 

As we make clear in Section 3.3, there is an intimate connection between the 
proximal-projection method and the well-known proximal point method with Breg- 
man distances - see p.23p . The proximal point method with Bregman distances 
considered in this paper is itself a generalization of the classical proximal point 
algorithm developed since the ninety-fifties by Krasnoselskii JEJ , Moreau [55] , [56] , 
[57] . Yosida [73], Martinet [53], [54] and Rockafellar |M|, [67| among others (see 
[50j for a survey of the literature concerning the classical proximal point method). 
It emerged from the works of Erlander 43 , Eriksson [42 , Eggermont 41j and Eck- 
stein [40 who studied various instances of the algorithm in M". Its convergence 
analysis in Banach spaces which are not necessarily Hilbertian was initiated in [30] , 
[3T] and [47] (see [32] and [30] for related references on this topic). Lemma 3.5 
shows that, in our setting, the proximal point method with Bregman distances is a 
particular instance of the proximal-projection method. 

Computing projections, metric or proximal, onto a closed convex set with com- 
plicated geometry is, in itself, a challenging problem. It requires (see ^^) solving 
convex nonlinear programming problems with convex constraints. Specific tech- 
niques for finding proximal projections are presented in [3, [21j and |33j . It is 
obvious from these works that it is much easier to find proximal projections onto 
sets with simple geometry like, for instance, hyperplanes, half spaces or finite in- 
tersections of such sets. These facts naturally led to the question of whether it 
is possible to replace in the process of computing iterates of metric or proximal 
projections algorithms the constraint set C by some approximations Ck of it whose 
geometry is simple enough to allow relatively easy calculation of the required metric 
or proximal projections at each iterative step k. That this approach is sound is quite 
clear from the works of Mosco [58] and its subsequent developments due to Attouch 
|15j . Aubin and Frankowska ^16j, Dontchev and Zolezzi 03, and from the studies 
of Liskovets [51], [52]. Its main difficulty in the case of the proximal-projection 
method is that the approximations Ck one uses should converge to C in a manner 
that ensures stable convergence of the algorithm to solutions of the problem. For 
some variants of the proximal-projection method, types of convergence of the sets 
Ck to C which are sufficiently good for this purpose are presented in [4], [7], [9], 
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[12j . [13] . They mostly are relaxed forms of Hausdorff metric convergence. It was 
shown in [8] that, in some circumstances, fast Mosco convergence of the sets Ck 
to C (see [5", Definition 2.1]), a form of convergence significantly less demanding 
than Hausdorff convergence, is sufficient to make the proximal-projection method 
p.7p applied to variational inequalities converge. As we show below, these con- 
vergence requirements in the case of the proximal-projection method (|1.7p can be 
further weakened. In fact, in our convergence theorems for the proximal-projection 
method we only require weak Mosco convergence of the sets Ck to C (see Defini- 
tion 4.3) and this is significantly demanding than Hausdorfl: metric or fast Mosco 
convergence. 

The purpose of this work is to find general conditions which guarantee that 
the proximal-projection method converges weakly or strongly to solutions of the 
Problem 1.1. Observe that there is no apparent connection between the data of 
Problem 1.1 and the function / involved in the definition of the proximal-projection 
method. Our main question is how the function / should be chosen in order to 
ensure (weak or strong) convergence of the proximal-projection method to solutions 
of Problem 1.1 without excessively conditioning the problem data. The function 
/ is a parameter of the proximal-projection method whose appropriate choice, as 
we show below, can make the procedure converge to solutions of Problem 1.1 even 
if the problem data are quite "bad" in the sense that they do not have some, 
usually difficult to verify in practice properties like maximal monotonicity, strict 
monotonicity, various forms of nonexpansivity, continuity or closedness properties 
of some kind or another. Theorem 4.1 and Theorem 5.1 are our responses to the 
question posed above. 

Theorem 4.1 shows that for guaranteeing that the proximal-projection method 
produces weak approximations of solutions for Problem 1.1 it is sufficient to chose 
a function / which, besides the conditions (|1.6p and (|1.8p which are meant to 
make the procedure consistent with the problem data, should satisfy some require- 
ments which, most of them, are common features of the powers of the norm ||-||^ 
with p > 1 in uniformly convex and smooth Banach spaces. The only somehow 
outstanding condition which we require for / is that it should be such that the 
operator A involved in the problem be Df-coercive on C or, if the set C is approx- 
imated by sets Ck, then D^-coercivity of A should happen on the union of those 
sets. Z?j-coercivity, a notion introduced in this paper (see Definition 3.2), is a 
generalization of the notion of firm nonexpansivity for operators in a Hilbert space 
(see (j4.34p ). Although in Hilbert spaces provided with the function / = ^ ||-|| this 
notion coincides with the notion of firm nonexpansivity and, also, with the notion 
of _D/-firmness introduced in [20] (see Definition 3.3 below), outside this particular 
setting the notions of Z?/-coercivity and _D/-firmness complement each other (cf. 
Section 3.2). In Section 4.3 we present several corollaries of Theorem 4.1 and exam- 
ples which clearly show that fitting a function / to the specific data of Problem 1.1 
may came naturally in many situations. If AT is a Hilbert space and ^ is a firmly 
nonexpansive operator, then the natural choice is / = ^\\-\\ (see Corollary 4.1). In 
this case the convergence of the proximal-projection method happens to be strong 
if the set of solutions of the problem has nonempty interior. If in Problem 1.1 we 
have C = A, then the problem is equivalent to that of finding a zero for the opera- 
tor V/ — V/ o Af, where Af is the ^/-resolvent of A (a notion introduced in 20J - 
see also p.lQp 'l and application of the proximal-projection method with Cfe = A to 
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V/ — V/o A/ is no more and no less than the proximal point method with Bregman 
distances mentioned above. The operator V/ — V/ o Aj is I?/-coercive whenever 
the operator A is monotone (cf. Lemma 3.4). This leads us to the application 
of Theorem 4.1 to the operator V/ — V/ o Af which is Corollary 4.2. It shows 
that the proximal point algorithm with Bregman distances converges subsequen- 
tially weakly (and when A has a single zero, sequentially weakly) for a large class 
of functions /, whenever A is monotone and provided that its graph is sequentially 
weakly closed (as happens, for instance, when Graphs is convex and closed in 
X X X*). It seems to us that this is the first time when weak convergence of the 
proximal point algorithm with Bregman distances is proved without requiring max- 
imal monotonicity of A. The proximal-projection method is also a tool for solving 
monotone variational inequalities via their Tikhonov-Browder regularization. This 
is shown by Corollary 4.3, another consequence of Theorem 4.1. Corollary 4.3 also 
asks for the monotone operator B : X -^ 2^ involved in the variational inequality 
to be such that V/ — V/ o Proj^ o [(1 — a) V/ — B] is I?/-coercive (or, equivalently, 
such that Proj^o [(1 — a)V/ — B] to be Df-Rrm). This happens in many situations 
of practical interest. Several such situations are described in the Examples 4.1 and 
4.2. 

A careful analysis of the proof of Theorem 4.1 reveals the fact that the proximal- 
projection method (|1.7p is a procedure of approximating fixed points for the opera- 
tor Proj^o(V/ — A) by iterating the operator. A customary tool of proving conver- 
gence of such algorithms in Hilbert spaces is the already classical Opial Lemma f59l 
Lemma 2]. Unfortunately, this result can not be extrapolated into a nonhilbertian 
setting in its original form. Our Theorem 4.1 is based on Proposition 4.1, a gen- 
eralization of the Opial Lemma which works in reflexive Banach spaces and which 
is of interest by itself. If the Banach space X has finite dimension. Proposition 4.1 
can be substantially improved - see Proposition 5.1. Thus, in spaces with finite 
dimension we can also improve Theorem 4.1 by dropping some of the requirements 
made on the problem data. This is, in fact, our Theorem 5.1 which guarantees con- 
vergence of the proximal-projection method with less demanding conditions than 
closedness of the graph of A. Accordingly, in finite dimensional spaces the con- 
clusions of Corollaries 4.2 and 4.3 can be reached at lesser cost for the operators 
involved in them as shown by Corollaries 5.1 and 5.2, respectively. 

This paper continues and develops a series of concepts, methods and techniques 
initiated in [T], [TS], [IHj, [SI], [3S] and [17] • In Sections 2 and 3 we present in 
a unified approach the notions, notations and preliminary results on which our 
convergence analysis of the proximal-projection method is based. It should be 
noted that, in Section 2, some of the notions and results are presented in a more 
general setting than strictly needed in the subsequent parts of the material. This is 
done so because we hope to use the framework created in the current paper as a base 
for a forthcoming study of methods of solving nonclassical variational inequalities 
which are only tangentially approached here. 

2. Proximal Mappings, Relative Projections and Variational 

Inequalities 

In this section we present the notions, notations and results concerning proxi- 
mal projections and variational inequalities which are essential for the convergence 
analysis of the proximal-projection method done in the sequels. 
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2.1. Proximal mappings and relative projections. All over this paper we 
denote by T^ the set of proper, lower semicontinuous, convex functions ^p : X ^ 
(—00, +00] which satisfy the conditions that 

(2.1) dom ip n int doni / 7^ 0, 
and 

(2.2) (/?f := inf \}p(x) : x € dovuLp n dom/} > — cxj. 

With every if e J-f we associate the function Env;^ : X* -^ {~oo, +cx)] given by 

(2.3) Env^(0 = miMx) + Wf{^,x) : x e X}. 

This is a natural generalization of the notion of Moreau envelope function (see |68[ 
Definition 1.22]). By (|2.ip and (|2.2p it results that the function Env;^ is proper 
and domEnv;^ = dom/*. Another generalization of the Moreau envelope function 
is the function env^ := Env^ o V/ introduced and studied in [22]. Using Fenchel's 
duality theorem, it is easy to deduce that if t/s € !Ff, then 

Env^(o = no -{v+ IT io = no " ('/'*□/*) (0- 

where ip*Df* denotes the infimal convolution of tp* and /*. 

The next result shows a way of generalizing the notion of Moreau proximal 
mapping (in the sense given to this term in [68]) whose study was initiated in 
[55] . [56] . [57] and further developed in [66], [67]. As we will make clear below, 
the generalization we propose here slightly differs from the notion of ^/-proximal 
mapping introduced and studied in ^20] . In fact, most of the next lemma can be also 
deduced from [20l Propositions 3.22 and 3.23] due to the equality (|2.7p estabhshed 
below. We prefer to present it here with a direct proof for sake of completeness. 

Lemma 2.1: Suppose that (p £ J^f. For any (, € int dom/* there exists a 
unique global minimizer, denoted Prox^(^), of the function (p{-) + W/(^, •). The 
vector Prox;^(^) is contained in dom 9(^ n int dom / and we have 

(2.4) Prox^(0 = [a (^ + f)]-' iO = {dip + V/)-i (0- 

Proof. Let £, G int dom/*. Then Env^(^) is finite and the function / — (^, •) is 
coercive (see [M] Theorem 7A] or [THl Fact 3.1]), that is, its sublevel sets 

lev{.(a) ■.= {x&X : fix) < a} , 

are bounded for all a g R. Consequently, the function Wf{£,, •) is coercive too. Let 
ja;'"'}, j^ be a sequence contained in domipD dom/ and such that 

hm [f{x'') + Wfitx'')] -Env^(0. 

The sequence {'/?(a;'^) + Vy/(^,x'^)}^ being convergent is also bounded. So, for 
some real number M > we have 

Wfi^, x'') < M - 99(0;''') < M - ipj, 

showing that the sequence {x'^} i.pm is contained in the sublevel set lev^(Af — 
LPf) of the function i/j := Wf {£,,■). By the coercivity of Wf{^,-) it follows that 
the sequence {2:''};.^^} is bounded. Since the space X is reflexive, it results that 
{x'^}^ has a weakly convergent subsequence {2;*'°}^.^^^ ■ Let x = w-limfe^oo a;'*" . 
The functions / and (p are sequentially weakly lower semicontinuous because they 
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are lower semicontinuous and convex. Hence, ip{-) + Wf {£,,■) is also sequentially 
weakly lower semicontinuous and, thus, we have 

ip{x) + Wfi^x) < liminf [ifix^") + Wfi^^x'")] = Env;^(C) < +00. 

This implies that x G domt/? n dom/ and that x is a minimizer oi ip{-) + Wf{^, •). 
Suppose that y is any minimizer of (/?(•) + Wf{^, •). Then y is also a minimizer of 
(f + f-^. Therefore, we have that e (9((^ + /-^)(y), that is, ^ € d{ip + f){y). The 
function / is continuous on intdom/ (as being convex and lower semicontinuous). 
This and (EH) imply (see [651) that d{ip + f){y) = d(p{y) + V/(y). Hence, 

(2.5) C e 5(^(2/) + V/(y). 
Since domV/ ~ intdom/ (see ()1.2p ). this implies that 

y G dom dip D dom V/ = dom dip D int dom /. 

Hence, all minimizers of ip{-) + W/(^, •) are contained in domdip D intdom/ C 
dom (p n int dom /. The Legendre function / is strictly convex on the convex subsets 
of dome}/ and, in particular, on the convex set doTatp n intdom/ = doimp n 
dom 9/. Thus, 'p{-) + Wf{S,, ■) is strictly convex on this set too. Consequently, there 
is at most one minimizer oi ip{-) + Wf ($, •) on the convex set dom cp (1 int dom / and 
this proves that the minimizer x whose existence was established above is unique. 
Formula (|2.4p follows from (|2.5p when y = x. D 

Lemma 2.1 ensures well definedness of the function 

(2.6) Prox;^ : intdom/* —^ dom9</jn intdom/, ^ -^ Prox^(^) 

when ip G J-^f. We call this function the proximal mapping relative to f associated to 
(p. Well definedness of the proximal mappings relative to the Legendre function / 
can also be deduced from [211, Theorem 3.18] where well definedness of the resolvent 

(2.7) prox;^ := Prox;^ o V/ 

was established. In more particular circumstances for / and ip, well definedness of 
Prox^ was proved in [1], [31], and in 36]. 



Let _E be a closed convex subset of X satisfying 

(2.8) ETlintdom/ 7^ 0. 

Then the indicator function of the set E, that is, the function le '■ X —^ (— oo, +oo] 
defined by le{x) =Q iix & E, and iE{x) = +oo, otherwise, is contained in Tf. The 
operator Prox^ is called projection onto E relative to / (cf. [36] ) and is denoted 
Proj;^ in that follows. According to Lemma 2.1, we have 

Proj{; - [Ne + V/)-i 
where Ne denotes the normal cone operator associated to the set E. The operator 

(2.9) proj;^ := proxf^ = Proj{ o V/ 

is exactly the Bregman projection onto E relative to f whose importance in convex 
optimization was first emphasized in [23] . To see that it is sufficient to recall that 
the Bregman distance Df : X x intdom/ — > (— cx),+cx)] is the function defined by 

(2.10) Df{y,x) := f{y) - f{x) - {Wf{x),y-x) = Wf{Wf{x),y), 
with dom 13/ = (dom/) x (intdom/) . 
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2.2. Variational inequalities. There is an intimate connection between proxi- 
mal mappings and variational inequalities. It is based on the following result which 
extends the variational characterization of Proj;^ given in |36j to a variational char- 
acterization of Prox^,. 

Lemma 2.2. Suppose that ip G Tf and ^ G intdom/*. If x d dom 9(p n 
int dom / then the following conditions are equivalent: 
(a) X =Prox/(C); 
(6) X is a solution of the variational inequality 

(2.11) {^-\/f{x),y-x) <ip{y)^ip{x), Vy G dom (p n dom /. 
(c) X is a solution of the variational inequality 

Wf{^,x) + Wf{Vf{x),y) - Wfi^, y) < ^(y) - ^(x), Vy edoimpn dom/. 

Proof: |(q)=^(o)| Suppose that x =Prox^(^). Take y G domiy9 n dom/ and t G 
(0, 1). Then (1 — t)x + ty E dom ip n dom / and we have 

V?(f) + Wf{^, x) < ^{{1 - t)x + ty) + Wf{^, (1 - t)x + ty), 

that is 

(2.12) ^((1 - t)x + ty) - ^{x) > Wf{^, x) - Wf{^, (1 - t)x + ty). 

Since x G intdom/, there exists to G (0,1) such that for any t G (0, to) we have 
that (1 — t)x + ty E int dom /. Consequently, for any t G (0, to) the function Wf {£,,■) 
is differentiable at (1 — t)x + ty. Clearly, we also have 

yWf{t-){u) = Vf{u)-^, Vu G intdom/. 

Therefore, by the convexity of W/(^, •) and (|2.12|) . we deduce 

t-^ [^{{1 ~ t)x + ty) ~ p{x)] > {\/Wf{t-){{l-t)x + ty),x-y} 

- (V/((l-t)£ + ty)-e,£-y) 

for any t G (0, to). The function V/(-) is norm to weak continuous (see, for instance, 
[5T1 Proposition 2.8]). Hence, letting t ^- O"*" in the last inequality we get 

ip°{x;y-x) > (V/(f) -C,x-y) , 

where ip° stands for the right-hand side derivative of ip, that is, 

Of n ,• ipix + sd) - ip{x) 
(p (x;a) ^ hm . 

s^0+ s 

Taking into account that 

(p{y) - (p{x) > (p°{x;y-x) 

we obtain ()2.11|) . 

\(b)^{a)\ Suppose that 

(C ^ "^ f{x),y — x) < p{y) — p{x), \fy G dom</7 n dom/. 

Observe that 

(2.13) VW^/(e,-)=V/(-)-e 
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Then, by the convexity of Wf{^, •) and ()2.13p . for any y e dom ip n dom/ we have 
that 

= (v/(£)-e,y-£) 

> ip{x) - ip{y). 

This shows that x =Prox^(f). The equivalence |(o)^(c)| results immediately by 
observing that 

Wf{^,x) + Wf{Wf{x),y)^Wf{^,y) = {^-Wf{x),y-x), 

whenever y e dom/, ^ G intdom/* and x G intdom/. D 

A consequence of Lemma 2.2 is the following generalization of the variational 
characterization of the Bregman projections originally given in [5]. 

Corollary 2.1. Let x € int dom/ and let E be a nonempty, closed and convex 
set such that E n int dom/ ^0. If x £ E, then the following conditions are 
equivalent: 

(i) The vector x is the Bregman projection of x onto E with respect to /; 

(ii) The vector x is the unique solution of the variational inequality 

(V/(x)~V/(z),z-y)>0, yyeE; 

(Hi) The vector x is the unique solution of the variational inequality 

Df{y,z) + Df{z,x)<Df{y,x), ^y e K. 

Now we are in position to establish the connection between the proximal map- 
pings and a class of variational inequalities. It extends similar results known to 
hold in less general settings (see, for instance, [l] and [45l Proposition 1.5.8]). The 
variational inequality we consider here is 

(2.14) Find x E intdom/ such that 

3^ e Bx : [{£„y-x) > (p(x) - (p{y), Vy e dom/] , 
where ip £ Tf and B : X ^ 2^ is an operator which satisfies the condition 

(2.15) 7^ domS n dom</3 n intdom/ and ran (V/ — B) C intdom/*. 
Condition (I2.15|) guarantees that the operator 

Prox^ (V/ - B) :- Prox^ o (V/ - B) 

is well defined. Therefore, the following statement makes sense. 

Lemma 2.3. Let if G ^f and x £ domdcpO intdom/. Suppose that B : X ^ 
2^ is an operator which satisfies (|2.15p . Then x is a solution of the variational 
inequality (|2.14p if and only if it is a fixed point of the operator Prox^ (V/ — B) . 

Proof: Note that x is a solution of (|2.14p if and only if there exists £, G Bx such 
that 

(2.16) ((V/(x)-- 0-V/(x),y-:E) <(p{y)-(f{x), Vyedom(^ndom/. 

According to Lemma 2.2, this is equivalent to a; = Prox^ (V/(a;) — £_) for some 
^ S Bx which, in turn, is equivalent to x € Prox^ (V/(a;) — Bx) , i.e., to the 
condition that x is a fixed point of Prox;;, {Vf — B) . D 
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Let B : X -^ 2^ be an operator and suppose that the closed convex subset C 
of X satisfies 

(2.17) 7^ dom S n C n int dom / and ran(V/-S) C intdom/*. 

Note that ii ip := ic, then the variational inequality (j2.14p is exactly a classical 
variational inequality 

(2.18) Find a:; G C n int dom / such that 

3^eBx:[{^,y-x)>0, Vy G Cndom/]. 

Applying Lemma 2.3 and (jl.3p in this case, we re-discover the following known 
resuh (cf. [2J): 

Lemma 2.4. Suppose that the condition ()2.17p is satisfied and that x (£ C D 
int dom /. Then the following statements are equivalent: 

(a) The vector x is a solution of the classical variational inequality (j2.18p ; 

(6) The vector x is a fixed point of the operator Proj^ (V/ — B); 

(c) The vector x is a zero for the operator V[B; C; f] : X -^ 2^ given by 

(2.19) V[B; C- f] := Vf~Vfo Proj^ (V/ - B) . 



Lemma 2.3 and Lemma 2.4 provided the initial motivation for this research. 
Observe that Lemma 2.4 reduces the problem of finding a solution for the classical 
variational inequality ()2.18|) to the problem of finding a fixed point for the operator 
Proj^ (V/ ~ B). It is well known that, in many instances, by iterating an operator 
starting from initial points located in its definition domain, one produces sequences 
which converge to fixed points of the operator. This suggests that for solving the 
classical variational inequality ()2.18p we would have to produce sequences defined 
by the iterative rule 

(2.20) x^+'^ e Proj^(V/(x'=) - Bx^) 

in hope that such sequences will converge to fixed points of Proj^ (V/ — B) . Note 
that any fixed point of the operator Proj^ (V/ — S) is a zero for the operator 
V[B\ C; f] given by (|2.19p and conversely. According to (|2.9p . we have that 

(2.21) Proj^(V/-y[B;C;/]) = Proj^ o V/ o Proj^ o (V/ - S) 

= proj^oproj^o (V/)"^o(V/-B) 

= proj^ o Vr o (V/ - B) 

= Proj^ (V/ - B) . 

Thus, we are naturally led to the question of whether, and in which conditions, the 
sequences generated according to the rule x'^'^^ G Proj^(V/(a;'^) — V[B;C; f]x''), 
which are the same (see (|2.2ip ) as the sequences generated according to rule (|2.20p . 
converge to zeros of the operator V[B; C; /]. This is, in fact, a particular instance 
of the more general question of whether, and in which conditions, the proximal- 
projection method (|1.7p approximates zeros of a given operator V[B;C;f], provided 
that such zeros exist. In the sequels we present answers to this question. It is in- 
teresting to observe that by focusing in our convergence analysis on conditions 
concerning the operator V[B; C; /] instead of B we do not mean that computing 



A METHOD FOR FINDING ZEROS OF SET-VALUED OPERATORS 11 

values of Proj^ (V/ — V[B; C; /]) is easier than computing values of the same op- 
erator via the formula Projp (V/ — B) . However, from a theoretical point of view, 
the operator V[B; C; f] associated to B via formula (|2.19p may happen to be better 
conditioned than B for a convergence analysis of the proximal-projection method. 
This aspect can be clearly seen after a careful dissection of the considerations which 
lead to our main convergence results presented in this paper. It should be taken 
into account that the operator B may have not zeros in C even if the operator 
V[B; C; /], associated to B by (|2.19p . has. For example, take X = R, f{x) — ^x^, 
C = [1, 2] and Bx — x for all x € X. Then V[B\ C; f]x = x — 1 vanishes at x = 1, 
but B does not have any zero in C. 

3. £)j-NONEXPANSIVITY POLES, Dj-COERCIVITY AND Dj-FIRMNESS OF 
OPERATORS IN BANACH SPACES 

In this section we introduce the notion of I?/-coercivity for operators from X to 
2^^ . We clarify how this notion is related with the notions of D/-nonexpansivity 
pole introduced in [32] and of D/-firm operator introduced in [2D]- Using these 
relations we show that the proximal point method with Bregman distances can 
be seen as a particular instance of the proximal-projection method applied to a 
-D/-coercive operator. 

3.1. Dy-nonexpansivity poles. In that follows, to the function / described in 
Section 1 and to any operator A : X —^ 2^ we associate the operator A^ : X ^ 2^ 
given by 

(3.1) A^:=Vro(V/-A). 

We call this operator the D f-antiresolvent of A. Observe that 

(3.2) dom A-^ C dom A n int dom / and ranA-^ C intdom/, 

and that, if x e intdom/, then 0* G Ax if and only if a; G FixA-^, where FixA-^ 
denotes the set of fixed points of A^ . Therefore, the (possibly empty) set of solutions 
of the Problem 1.1 situated in intdom/, denoted Sf{A,C), is exactly 

(3.3) Sf{A,C) = CC^VlliAf. 

We are going to prove that, for operators A which are Z?/-coercive, the set Sf{A, C) 
is exactly the set of Dy-nonexpansivity poles of A^ over the set C and this fact will 
be later used in our convergence analysis of the proximal-projection method. To 
this end, recall the following notion. 

Definition 3.1 (cf. [32]) Let T : X ^ 2"^ be an operator and let F be a subset 
of X such that 

(3.4) ^ T{Y n int dom /) C int dom /. 

The vector z ^ X ia called a D f-nonexpansivity pole of T over Y if the following 
conditions are satisfied: 

(3.5) zeFn intdom/, 

(3.6) (a; e r n int dom / and u eTx) ^ {V f{x) -V f{u),z ~~ u) <{). 
We denote by Nexpy T the set of _D/-nonexpansivity poles of T over Y. 



12 DAN BUTNARIU AND GABOR KASSAY 

Operators having Z?/-nonexpansivity poles were termed totally nonexpansive 
operators in |32j . Operators T such that ran T C dom T = int dom / and hav- 
ing Nexp/^^ j^jj^ . T D FixT were called B -class operators in [12] and [2D|. ;B-class 

operators necessarily have Nexp/^^. ^^^^ . T = FixT (cf. [20l Proposition 3.3]). How- 
ever, not every operator having I?j-nonexpansivity poles over some subset Y of 
X is S-class. For example, the operator Tx — {a;^} when X = M, / = ^ |-| and 

Y = int dom / = M has z e Nexp/^j ^^^ , T if and only if 

x{l-x){z-x^) <Q, Va;eR, 

and this last inequality can not hold because, irrespective of z, we have lima;_,oo a;(l— 
x){z — x"^) = cx). Hence, Nexpfjj^ ^^^^^ ^ T = 0. In spite of that, FixT = {0, 1} and, if 

Y = [0, 1], then dSH) holds for z = only, i.e., Nexp^^^^j T ^ {Q} ^ FixT. 

The following lemma summarizes several properties of operators having nonex- 
pansivity poles which are used in this work. 

Lemma 3.1. Let the operator T : X —^ 2^ and the set Y C X be such that 
condition (|3.4|) holds. Then the following statements are true: 

(a) The (possibly empty) set Nexp^ T is convex and closed when Y C int dom / 
is convex and closed; 

(b) A vector z (z Y int dom/ is a D f -nonexpansivity pole of T over Y if and 
only if 

(3.7) (a: e F n int dom / and u E Tx) ^ Df{z,u) + Df{u,x) < Df{z,x). 

(c) Nexpy T CFixT and T is single-valued at any Df -nonexpansivity pole. 

Proof. Statement (a) results from the fact that the function z -^ 
(V/(a;) — Wf{u),z — u) is linear and continuous of z. Statement (6) follows from 
((3:6)) and (|2J0)) . Now, by taking in dSJ]) x = z e Nexp(, T one gets 

(3.8) Df{z, u) = Df{u, z) = 0, Vu e Tz. 

By p.4|) . if w £ Tz, then u E int dom/. The function / being essentially strictly 
convex is strictly convex on int dom /. Hence, the equalities in p.Sp can not hold 
unless u = z (cf. [32l Proposition 1.1.4]). In other words, we have the following 
implication 

(3.9) zG Nexp(,T^Tz = {z}. 

f . 

It shows that Nexpy T C Fix T and that T is single- valued at D^-nonexpansivity 

poles. D 

3.2. Dy-coercivity and D^-firmness. The notion of Z?/-coercive operator, in- 
troduced in this section, and the notion of _Dy-firm operator, originally introduced 
in [201 Definition 3.4], are generalizations of the notion of firmly nonexpansive op- 
erator in a Hilbert space. Recall (cf. [46l pp. 41-42]) that if X is a Hilbert space 
(which we always identify with its dual X*), then an operator A: X -^ X is firmly 
nonexpansive on a subset 1" of X if and only if 

(3.10) {Ax- Ay,x~y) > \\Ax-Ayf, ^x^y E Y. 

It can be easily seen from the definitions given below that if X is a Hilbert space 
and if /= 5 II' II , then the operator A is T)/-coercive on its domain if and only if it 
is firmly nonexpansive on its domain and this happens if and only if A is D/-flrm. 
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Returning to the general context in which X and / are as described in Section 
1, we introduce the following notion. 

Definition 3.2. Let F be a subset of the space X. The operator A : X ^ 2-^ 
is called Df- coercive on the set Y if 

(3.11) rn (doniA)n (intdom/) 7^0 

and 
(3.12) 

x,y GY niiit dom / 



C e ^. and , . A, r - ( e - ., V/* (V/(.) - - V/* (V/(.) - ,)) > 0. 

Operators satisfying a somewhat less restrictive condition than (|3.12p were stud- 
ied in [36[ Section 5] under the name of inverse-monotone operators relative to /. In 
general, an operator A (even in a Hilbert space provided that / is not the function 
/ = I ll'll ) does not have to satisfy p.lOp in order to be D/-coercive on Y. For in- 
stance, if the function / has a minimizer in int dom/ (i.e., if the equation V/(a;) = 
has a solution), and if a € (0,1), then the operator A = aV/ is Z3j-coercive on 
Y — dom A = int dom/ without necessarily satisfying p.lOp on Y ~ int dom/. 
Indeed, in this case, ii x £ int dom/ then {V/(a;),0*} C ranV/ = int dom/* and, 
due to the convexity of int dom /* , we have that 

(1 - a)Vf{x) = (1 - a)Vf{x) + aO* G int dom/* = dom V/*, 

and, consequently, the operator 

Afx = Wr{{l-a)Wf{x)) 

has doniA-'^ = int dom/. Ii x,y E int dom/, and ii (3 = 1 — a, then, by the 
monotonicity of V/*, we deduce that 

( Ax - Ay, V/* (V/(x) - Ax) - V/* (V/(y) - Ay)) = 
a/3-1 (^vj(2,) _ ^v/(2/), V/*(/3V/(a;)) - V/*(/3V/(y))) > 0, 
i.e., the operator A = aV/ is £)j-coercive on int dom/. However, the operator A = 
aVf does not have to be firmly nonexpansive on its domain even if X is a Hilbert 
space. For example, take in the considerations above X — L^[0, 1], / = i ||-|| and 
a G (1/4, 1). Then V/(a;) = ||a;|| x and an easy verification shows that (|3.10p does 
not hold for any x,y £ X. For instance, p.lOp is violated when x = 2q~^I'^ and 
2/ = 0. 

We are going to show that there are strong connections between the D^-coercivity 
of the operator A and the D/-firmness of its D^-antiresolvent A^ . For this purpose 
we recall the following: 

Definition 3.3. (cf. [SO", Definition 3.4]) An operator T : X ^ 2^ \& called 
Df-firm if it satisfies the conditions 

(3.13) 7^domrUranrcintdom/ 
and 

(3.14) ueTx and veTy=> (V/(u) - V/(w), u ~ v) < (V/(a;) - V/(y), u - v) . 

We start with the following result which summarizes some basic properties of 
Z^Z-coercive operators. 

Lemma 3.2. Let A : X ^ 2^ be an operator and let Y be a subset of X 
which satisfies p. lip . The following statements are true: 



£ Ah' ( -^ Df{u,v) + Df{v,u) + Df{u,x) + Df{v,y) < Df{v,x) + Df{u,y); 



14 DAN BUTNARIU AND GABOR KASSAY 

(a) The operator A is Df -coercive on Y if and only it satisfies the following 
condition for any x,y £ Y D int dom / : 

u G A^x 
V e A^y 

(b) If A is Df -coercive on Y, then A is Df -coercive on any subset of Y which 
intersects dom A n int dom /; 

(c) The operator A is Df -coercive on its domain if and only if its D f -antiresolvent 
A^ is Df-firm; 

(d) // A is Df -coercive on its domain, then A^ is single valued and all its fixed 
points are D f-nonexpansivity poles on int dom/. 

Proof. Statements (a), (6) and (c) result from ([21]), dSSD, ((XTO)l . ((XTI)) and 
(|2.10p . Single valuedness of A^ in statement (d) is a consequence of (c) and of 
[20l Proposition 3.5(iii)]. Letting y = z £ Fix^-'' in the inequality of (a), and 
taking into account the single valuedness of A^ , one obtains that z satisfies p.7|) 
{0TT = Ah D 

Whenever the operator A involved in Problem 1.1 is Z?/-coercive on the domain 
C of the problem, we have 

(3.15) Sf{A, C) = C n Fix A^' = Nexp^ Ah 
This immediately follows from the next result. 

Lemma 3.3. The following statements are true: 

(a) If ze Nexp^ A^ , then Az = {0*}; 

(b) If the operator A is Df -coercive on C and z e int dom/ is a solution of 
Problem 1.1, then z G Nexp^ A-^. 

Proof. The statement (a) results from (jS.ip and Lemma 3.1(c). In order to 
prove (&), note that for any x £ C C\ int dom/ and y e A^ x we have 

Df{z, x) - Df{z, y) = Df{y, x) - {Wf{x) - Wf{y),z - y) , 

and y = Vf*{Vf{x)- C) for some ^ S Ax. Thus, Vf{y) = V/(a;)- ^ and 

(3.16) Df{z, x) - Df{z, y) = Df{y, x) + (0* - C, ^ " 2/> ■ 

If z is a solution of Problem 1.1, then z G C, 0* G Az and z = Wf*{Wf{z) - 0*). 
By consequence, 

(3.17) (0* - ^, z - y) - (0* - e, ^f*{^f{z) - 0*) - V/*(V/(a;) - 0) ■ 

Since A is D/-coercive on C, it results that the right-hand side of (|3.17[) is non- 
negative (see p.l2p ) for any x £ C (1 int dom/. Hence, by p.l6|) and (j3.17|) . the 
inequality in p.7p results and it shows that z e Ncxp^ A^ . D 

The class of operators which are Z?j-coercive contains some meaningful opera- 
tors. Among them are all operators A[T] : X ^ 2^ given by 

(3.18) A[T]=V/-V/or, 

where T : X — > 2^^ is a D/-firm operator. In particular, A\T\ is D/-coercive when 
T = Bf, where Bf : X ^ 2-^ is the Df -resolvent (cf ^20^) of a monotone operator 
B :X ^ 2^', i.e., 

(3.19) Bf.^iVf + Br'oVf. 
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According to [20l Proposition 3.8], the operator T = Bf satisfies the condition 
()3.13|) . These facts are summarized in the fohowing lemma. 

Lemma 3.4. Let T : X -^ 2^ be an operator which satisfies (|3.13p . Then the 
following statements are true: 

(a) dom A[r] = domT and its antiresolvent is A [T] = T; 

(b) The operator T is Df -firm if and only if the operator A[T] : X ^ 2^ 
defined by (j3.18p is Df -coercive on its domain. 

(c) If B : X ^ 2^ is a monotone operator with doniB n intdom/ 7^ 0, then 
Bf is Df-firm, single valued on its domain, 

(3.20) A[B/]-V/-V/o(V/ + B)-'oV/, 
and A[Bf] is D f-coercive on its domain. 

Proof. Statement (a) results from (|3.ip and (|3.18p . To prove [b) observe that 
for any x,y € domT, for any ^ € A[T]x and for any rj Cz A [T] y , we have 

(3.21) ^ = Vfix)-Vf{u) and ry - V/(y) - V/(i;) 
for some u G Tx and for some v G Ty. Therefore, 

(3.22) ( e - r,, V/*(V/(a;) - C) - V/*(V/(2/) - r;)) 

= {Vf{x) - Vf{y),u -v)^ (V/(u) ~ Vf{v),u ~ v) . 

If T is Df-finn, then the right hand side of p.22p is nonnegative and this implies 
that A [T] is _D/-coercive on domT. Conversely, suppose that A [T] is T)/-coercive 
on domT. If u G Tx and v G Ty, then the vectors f and 77 given by p.2ip satisfy 
(I3.22P and the left hand side of this equality is nonnegative. Hence, T is T'y-firm. 
This proves (6). In order to prove (c) recall that, since B is monotone, its resolvent, 
Bf, is necessarily Df-Hrin and single valued on its domain (cf. [201 Proposition 
3.8]) and, as noted above, it satisfies ()3.13p . Thus, according to (b), the operator 
A[i?/] is Z?/-coercive on its domain. D 

3.3. Connection between the proximal-projection method and the prox- 
imal point method. Lemma 3.4 helps establishing a connection between the 
proximal-projection method and the proximal point method (with Bregman dis- 
tances). The proximal point method we are referring to in this paper is the iterative 
procedure which, in our setting, can be described by 

(3.23) x°edomBf and x'^+i == B/(a;'=), V/c G N, 

where B : X ^r 2^ is a monotone operator with domi? n intdom/ 7^ 0. Its well 
definedness is guaranteed when 

(3.24) 7^ ranB/ C domS/. 

For ensuring that the inclusion in this condition holds it is sufficient to make 
sure that domi?^ = X. This implicitly happens when one considers the classi- 
cal proximal point method (see [6^) where X is a Hilbert space, / = \\\-\\ and 
B is presumed to be maximal monotone. Alternative conditions which imply that 
AoxQ. Bf = X when B is maximal monotone are presented in [30J in a more gen- 
eral setting. In particular, those conditions hold if A" is a uniformly convex and 
uniformly smooth Banach space, B is maximal monotone and / = 5 ||'|| . 

Maximal monotonicity of i? is a commonly used condition for ensuring well de- 
finedness of the proximal point method because, if B is maximal monotone and 
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/ = ^ IMI^ , then Vf + B and V/* are surjective and, thus, V/* o (V/ + B) is 
surjective too, that is, dom Bf — ran [V/* o (V/ + S)] = X. However, well de- 
finedness of the proximal point method can be sometimes ensured for operators 
B which are monotone without being maximal monotone. In such cases, it is in- 
teresting to know whether the proximal point method preserves the convergence 
properties which make it so useful in applications requiring finding zeros of max- 
imal monotone operators. Here is an example of a monotone operator which is 
not maximal and for which p.24p holds (in spite of the fact that V/ + B is not 
surjective). Take X ^ R, f = ^\f and let B : X ^ 2^' be given by Bx = {0} if 
a; < 0, and Bx = if x > 0. The operator B is monotone, but it is not maximal 
monotone since the operator defined by B'x = {0} for all x G X is a proper mono- 
tone extension of B. Obviously, in this case V/ is the identity, (V/ + B)x — {x} if 
X <Q, {Vf + B)x = if X > and, therefore, V/ -f i? is not surjective. However, 
Bf = {Wf + B)~ = Wf + B and, hence, ran_B^ — domBj — (— oo,0] showing 
that (|3.24p is satisfied, that is, the proximal point algorithm is well defined. 

The next result establishes the connection between the proximal-projection method 
and the proximal point method. It requires that dom Bf should be convex and 
closed. This necessarily happens if V/ + i? is surjective and dom/ = X. However, 
dom_B^ may happen to be convex and closed even if V/ -|- -B is not surjective as 
one can see from the example above. Other instances in which domBf is convex 
and closed are described in the remarks preceding Corollary 4.2 as well as in the 
body of that corollary. 

Lemma 3.5. Let B : X ^ 2^ be a monotone operator such that doraBf 
is convex and closed and suppose that (j3.24p is satisfied. Then the proximal point 
method p.23p is exactly the proximal-projection method applied to the Df -coercive 
operator A[Bf] with Ck = X for all fc S N. 

Proof. Observe that, by (|3.20p . we have 



;/ fY7 f _ 4rn.n — v>^,,J 



P'-OJ^omB,(V/-A[S/]) - PrOJLnS, V/ O (V/ + i?)"^ O V/ 



\-l 



ProJdo^B^ 



(V/ + B) 'oV/ 



Bf, 



where the last equality results holds because p.24p is satisfied. This shows that the 
proximal point method and the proximal-projection method are overlapping when 
one takes A = A[Bf] and Ck ^ X for all fc e N in (fL7| . D 

4. Convergence Analysis of the Proximal-Projection Method 

In this section we present a convergence theorem for the proximal-projection 
method in reflexive Banach spaces. Our convergence analysis is based on a gener- 
alization of Lemma 5.7 in >36\ which, in turn, is a generalization of a result known 
as Opial's Lemma [SSI Lemma 2] . All over this section we assume that the function 
/ and the Banach space X are as described in Section 1. 

4.1. A generalization of Opial's Lemma. Opial's Lemma says that if X is a 
Hilbert space and ii T : Y ^ X is a nonexpansive mapping on the nonempty 
closed convex subset Y of X, then for any sequence {z''} , ., C Y which is weakly 



is 



convergent and has limfe^oo \\Tz'' ~ -^'^11 ~ 0' ^^^ vector z = w — limfe_,oo z 
necessarily a fixed point of T. In [36|, Lemma 5.7] a similar result was shown to hold 
in Banach spaces which are not necessarily Hilbert spaces. Namely, it was proved 
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that the conclusion of Opial's Lemma still holds for operators T : X -^ X which 
are nonexpansive relative to / (in the sense given to this term in [37J), i.e., such 
that 

(4.1) Df{Tx,Ty)<Dfix,y), ^x,y eY, 

provided that dom/ = domV/ = X and that / it is not only Legendre, but it 
is also totally convex (see [32]) and bounded on bounded sets, while T satisfies 
limfc^oo%(^z^z^■)=0. 

Our current generalization of Opial's Lemma concerns set-valued operators T 
satisfying a somehow less stringent nonexpansivity condition than (j4.ip with re- 
spect to a function / subjected to weaker requirements than those involved in [36[ 
Lemma 5.7]. In the sequels we use the following notion which generalizes that of 
nonexpansive operator relative to /. 

Definition 4.1. The operator T : X — > 2^ is said to be Df -nonexpansive if it 
satisfies p.l3|) and for any x G dom T there exists u £ Tx such that 

(4.2) (Vy e domT) : [v G Ty^ Df{v,u) < Df{y,x)] . 

In that follows (see Theorem 4.1 below) we will be interested in operators whose 
antiresolvents are simultaneously D/-firm and I?/-nonexpansive. It should be noted 
that the notions of Dy-nonexpansivity and Z?/-firmness are not equivalent, although 
some operators may have both properties. If X is a Hilbert space provided with 
/ = i ||-|j , then it is obvious that any Dj-firm operator (i.e., any firmly nonexpan- 
sive operator) is 13/ -nonexpansive (i.e., nonexpansive). Even in this context, the 
converse implication does not generally hold. Take, for example, the case where 
the Hilbert space is X = R and Ax = 2x. Its antiresolvent is A^x = —x and 
it is 13/ -nonexpansive without being Z3/-firm. However, operators whose antire- 
solvents are simultaneously I3/-nonexpansive and I3/-firm are not specific to the 
setting of Hilbert spaces provided with / = | ||-|| . For instance, if X = M and 
f{x) = jx'^, then the antiresolvent A^ of the operator A = aS/ f with a £ (0, 1) is 
Z3/-nonexpansive and Dj-firm at the same time. It is ZPy-firm because, as shown in 
Section 3.2, A is D/-coercive and Lemma 3.2(c) applies. A^ is also I3/-nonexpansive 
because 

Df{Afy,Afx) = (1 - afDf{y,x) < Dfiy,x), Va:,y G X. 

One still may hope that (as happens in the particular situation noted above when 
X is a Hilbert space provided with / = 5 |H| ) a D/-firm operator is always Df- 
nonexpansive. The following example shows that this is not the case. 

Example 4.1. A Df-firm operator which is not D f -nonexpansive. Let X = K 
and let / : M ^- M be the Legendre function given by f{x) — |a;p/^. Take the 
continuous, single-valued operator T : M -^ K. defined by 

\^/x, ifa;>^. 
We first show that T is _D/-firm, that is, we verify that for every x,y E M. one has 

(4.3) ifiTix)) - nnyMTix) - Tiy)] < [f'{x) - f'iyWix) - Tiy)]. 

For symmetry reasons we can assume that x > y. The case x,y < 1/16 being 
trivial, we shall consider the following two cases. 
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Case 1. x,y > 1/16. In this case ()4.3|) reduces to [\/^ — ^][y/x — y/y] < 

{y^ ~ \/yY , which is equivalent with the obviously true inequality 1 < -y/ir + ^. 
Case 2. X > 1/16 and y < 1/16. In this case we distinguish two subcases: 
(i) y > 0. In this situation (|4.3p can be re- written as [\/x — yJ\/4\[-\/x — 1/4] < 

[^Jx — Y^][\/i' — 1/4], which, in turn, is equivalent with -y/ir + ^ < y/x + 1/2. 

This last inequality holds since for x > 1/16 and y < 1/16 we obviously have 

-^ + y/y < \px + 1/4 and it is easy to verify that 

(4.4) ^ + 1/4 < V^ + 1/2, Vx e [0, c»). 

(ii) y < 0- In this case, the inequality (|4.3p is equivalent to [\/x — \/l/4][y^ — 
1/4] < [y/x + y/^[y/x~l/4:], that is, to ^~l/2 < y/x + y/^. This last inequality 
is true because ^/x + ^J~^ > \fx and, by (|4.4p . we also have \fx < ^Jx + 1/4 < 
y/x + 1/2 for every a; > 0. 

These show that the operator T is ZJy-firm. Now we verify that T is not Df- 
nonexpansive, that is, that there exists two real numbers x and y such that 

(4.5) j{T{y)) - !{T{x)) - !\T{x)){T{y) - T{x)) > f{y) - f{x) - f{x){y - x) 

Taking x = 1/8 and y ~ 1/16, a simple computation shows that the left hand side 
of g31) equals 2^3/2 [23/4 ^ 31/2 _ 3] ~ 8.6895, while the right hand side equals 
2i5/4[2i/2_i]^ 5.5730. D 

Before proceeding with the presentation of our generalization of Opial's Lemma 
several observations concerning its hypothesis are in order. 

Remark 4.1. (a) The next result is a proper generalization of Lemma 5.7 in [36] 
which, in turn, is a proper generalization of Opial's Lemma. To see this, observe 
that when the operator T : X ^ X is Dy-nonexpansive, it satisfies ()4.7p below. 

ih) The hypothesis of the next result implicitly require that the space X should 
be reflexive. The fact is that a function / which is lower semicontinuous with 
int dom / 7^ and uniformly convex on bounded subsets of int dom / exists on a 
Banach space X only if X is reflexive (cf. [34j Corollary 4.3] in conjunction with 
[361 Theorem 2.10(u)]). 

With these facts in mind we now proceed with the presentation of the general- 
ization of Opial's Lemma. 

Proposition 4.1. Suppose that the function f is uniformly convex on bounded 
subsets of int dom /. Let T : X — > 2-^ be an operator satisfying p.l3p and suppose 
that V/ is bounded on bounded subsets of domTU ranT. // {z'^j^gN C domT is 
a sequence which converges weakly to a vector z G domT and if, for some sequence 
{"''IfceN satisfying 

(4.6) (V/c e N : u*' e Tz'') and lim Df{u'', z'') = 0, 

k — ^CXD 

there exists u G T{z) such that 

(4.7) liminfD/(M'=,M) < liminf ^/(z'', z), 

k — *oo k — s-oo 

then the vector z is a fixed point of T . 
Proof. For any x G int dom / one has 

Df{z\ x) - Df{z\ z) = Dfiz, x) + (V/(x) - V/(z), z-z''). 



A METHOD FOR FINDING ZEROS OF SET- VALUED OPERATORS 19 

This implies that 

(4.8) \immiDf{z'',x) > Df{z,x)+limMDf{z'',z), 

fc— *-oo A;— >oo 

because, since {^^} ^^^ converges weakly to z,we have 

hm(V/(x)-V/(z),z-z'=)=0. 

k — ^oo 

Since the function / is Legendre, it is strictly convex on int doni /. This implies 
that Df{z,x) > whenever x ^ z (cf. [32l Proposition 1.1.4]) and, consequently, 
by (l48l) we obtain 



(4.9) X ^ z ^ \imiiiiDf{z'',x) > liminf D/(2'', 2). 

k — ^00 k — ^00 

We claim that 

(4.10) limmfDf{u'',u) =\immi Df{z'',u). 

k — »oo A; — *oo 

To prove this claim, observe that 

(4.11) Df{u\u) = Df{z\u) + [f{u'^) - f(z^)] - {Vf{u),u'^ - z^). 

The sequence {^'^j^.^j^ is bounded as being weakly convergent. Since the func- 
tion / is uniformly convex on bounded subsets of intdom/, it is also sequentially 
consistent (cf. 36, Theorem 2.10]). Therefore, by (|4.6p . we deduce that 

(4.12) lim \\u'' - z''\\ =0. 



Hence, {w'^j^, pj is bounded and 



I fee 

(4.13) lim (V/(u),u''-z'=) =0. 

fc— »oo 

The convexity of / on int dom / implies 

(4.14) {'^f{u''),u^ ~z'')>f{u^)- f{z'')>{Vf{z^),u^ ~z^), VfceN. 

By hypothesis, V/ is bounded on bounded subsets of doniTU ranT. Therefore, 
the sequences {V f{z^y}k£n and {V f {u'^)} ke^ are bounded. Thus, by (|4.12p and 
(|4.14p . we deduce that 

hm [/(w'^) - /(z^-)] =. 0. 

k — >oo 

This, combined with (|4.1ip and (|4.13p . implies (|4.10p and the claim above is proved. 
Suppose by contradiction that z ^ Tz. Then the vector u e Tz whose existence 
is guaranteed by hypothesis has u ^ z and then, by (|4.9p . we deduce 

(4.15) liminfL>/(z'',u) > liminf D/(z'=, z). 

k — *oo k — *oo 

On the other hand, by (|4.7p and (|4.10p . we have that 

(4.16) liminf Z?/(z'',z) > lim inf L»/ (w'' , u) = liminf L»/(z'', w), 

k—^oG k—>-oo k—^oo 

which contradicts (|4.15[) . This completes the proof. D 
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4.2. A convergence theorem for the proximal-projection algorithm. At 

this stage we are in position to consider the question of convergence of the procedure 
(|1.7p towards solutions of Problem 1.1. For this purpose, we recall the following: 

Definition 4.2. (Cf. [58j ) (a) The weak upper limit of the sequence {Ek}i^^jq 
of subsets of X is the set denoted w-lhnk^ooEk and consisting of all a; e A such 
that there exists a subsequence {-Ej^l^gN of {-Efej^gN ^^'^ ^ sequence {x^} f^^^^ in 
A which converges weakly to x and has the property that x'' G Ei^ for each fc G N. 

(6) The operator A : X ^r 2^ is sequentially weakly- strongly closed if its graph 
is sequentially closed in A x A* provided with the {weak) x (siron(7)-topology, that 



4-17 fc ^ \ ,k , \^iGAv. 

V ^- V and ^ -^ ^ J 

Before proceedings towards the main result of this paper the following observa- 
tions may be of use. 

Remark 4.2. (a) The subset i? of A may happen not to be convex even if 
there exists a sequence of closed convex sets {-E/c}fcgN contained in A such that 
w-limfc^oo-Efe = E. Indeed, take A = R and 



^ fc+l' ^ ' k+l 

Ek ' 

I 1 , , 

fc+l' ' fe+1 



if k is even. 



if k is odd. 



It is easy to verify that the weak upper limit of this sequence of closed convex sets is 
the nonconvex set {—1, +1} . This fact explains why, in the next theorem, convexity 
of C can not be derived from the convexity of the sets Ck ■ 

(b) Any nonempty closed convex subset i5 of A is the weak upper limit of a 
sequence of half spaces (corresponding to support hyperplanes) containing it. 

(c) Among the operators which are sequentially weakly-strongly closed are all 
the maximal monotone operators (see, for instance, |60j). 

{d) An essential part of condition (6) of the theorem below is the requirement that 
the gradient V/ of the Legendre function / should be bounded on bounded subsets 
of intdom/. If the Legendre function / has the property that V/ is bounded 
on bounded subsets of intdom/, then dom/ = A. Indeed, since / is essentially 
smooth, it follows that intdom/ ^ and for any sequence {a^'^/i.pRj contained in 
int dom / and converging to a point of the boundary of int dom / has the property 
that Hmfe^oo ||^/(2;'^)||, = oo (cf. [TH Theorem 5.6]). Now, suppose that {a;'''}j,gp, 
is a convergent sequence contained in intdom/ and denote by x its limit. Then 
the sequence {V/ (a^*^)}/.;;]!^ is bounded because the sequence {x''},^^^ is bounded 
and V/ is bounded on bounded subsets of intdom/. We claim that x € intdom/. 
Assume by contradiction that x ^ int dom /. Then x belongs to the boundary of 
intdom/. Hence, lirufc^oo V/(a;'^) = oo and this contradicts the boundedness 
of {V/ (2;*') /i.pRj • Since intdom/ contains the limit of any convergent sequence of 
vectors contained in it, it follows that intdom/ is, simultaneously, a closed and 
open set. The space A, being a Banach space, it is necessarily arcways connected 
and, thus, a connected space (cf. [351 Theorem 10.3.2]). Consequently, A is the 
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only nonempty subset of X which is open and closed at the same time (cf. [381 
Theorem 10.1.8]), that is, intdom/ = X. 

The following theorem establishes the basic convergence properties of the proximal- 
projection method. It should be observed that, in view of Remark 4.2(d), the hy- 
pothesis of point (6) of the theorem implicitly requires that dom f = X. At point 
(ii) of the theorem sequential weak-weak continuity of V/ is mentioned as a suf- 
ficient condition for weak convergence (as opposed to subsequential convergence) 
of the proximal-projection method. This condition is obviously satisfied whenever 
the space X has finite dimension. It is also satisfied if X is a Hilbert space and 
/ = ^ ||-|| as well as in some nonhilbertian Banach spaces like P provided with 
/ = i ||.||P for any p G (1, +oo) - see [2H1 Proposition 8.2]. 

Theorem 4.1. Suppose that the function f is uniformly convex on bounded 
subsets of intdom/, V/* is bounded on bounded subsets of V/(domA) and that, 
in addition to (jl.6p and Assumption 1.1, the subsets C'k of X satisfy 



(4.18) C = w-limfc^ooCfc. 

If Problem 1.1 has at least one solution, if the operator A is Df -coercive on the set 
Q := Ufcgpj Cfc, and if at least one of the following two conditions is satisfied: 

(a) V/ is uniformly continuous on bounded subsets of intdom/ and A is se- 
quentially weakly-strongly closed; 

(b) V/ is bounded on bounded subsets of intdom/, A^ is D f -nonexpansive and 
C C dom A; 

then any sequence {x^^k&i generated by the proximal-projection method (|1.7p has 
the following properties: 

(i) It is bounded, has weak accumulation points and any such point is a solution 
of Problem 1.1; 

(a) If Problem 1.1 has unique solution or if \7 f is weakly-weakly sequentially 
continuous, then the sequence {x''}k£fi converges weakly and its weak limit is solu- 
tion to Problem 1.1; 

(Hi) If the Banach space X has finite dimension, then {x'^'j^gN converges in 
norm to a solution of Problem 1.1. 

Proof. Let z € C be a solution of Problem 1.1. Then, clearly z G dom A and, 
by (|1.6p . we deduce that z € intdom/. For each fc € N, let C'^ G Ax'' be such that 

(4.19) x^-+i = Proj^,^,ndomA(V/(x^-) - C')- 
Denote by 

(4.20) M^:=V/*(V/(x'=)-C''). 
Observe that 

(4.21) u'^eAfx'' and x'=+' = proj{,^^^ndomX, Vfc G N. 

By hypothesis, the operator A is Z?/-coercive on the set Q and z is a solution of 
Problem 1.1. Note that, since 0* € Az and z G C C Q, Lemma 3.3 applies with C 
replaced by Q. It implies that z sNexpg A^ . By Lemma 2.1 we have that 

(4.22) x*-' e CfcH dom An intdom/ C gn dom An intdom/, Vfc G N. 
First we prove the following: 

Claim 1: The sequence {2;'^}j,pR, is bounded. 



22 DAN BUTNARIU AND GABOR KASSAY 

In order to show this notice that, by applying Lemma 3.1(6) to z GNexp^ A^ we 
deduce 

(4.23) Dj{z,u^)^Df{u^,x^)<Dj{z,x^), Vfc G N. 
This imphes 

(4.24) Df{z,u'') <Df{z,x''), Vfc G N. 

By Assmnption 1.1, we have that z E C C Cfe, for aU fc e N. Thus, taking into 
account ()2.10|) . (I4.21|l and Lemma 2.2 apphed to ip — ic'k+indomA, we obtain that 

(4.25) Df{z,x''+^) + Df{x''+'^,u'')<Df{z,u''), Vfc G N. 
Combining (|i^ with ((i?^ yields 

(4.26) Df{z,x''+'^)<Df(z,x''), Vfc e N, 

showing that the nonnegative sequence {Df{z, a:;'°)}fcgN is nonincreasing and, there- 
fore, bounded. Let (3 be an upper bound of {Df{z,x'')}kefii- According to (II. 4|) . 
(f2l0l) and ([4261) . we deduce that 

n\7f{x''j) - (V/(a;'=), z) + /(z) = WfiWfix'^), z) = Df{z, x^) < (3, Vfc G N. 

This implies that the sequence {V/(a;'^)}fegN is contained in the sublevel set 
lev^ (/? — f{z)) of the function tp :— /* — (•, z) . Since z £ intdom/ = intdom (/*)*, 
and since the function /* is proper and lower semicontinuous, application of the 
Moreau-Rockafellar Theorem ([Ml Theorem 7(A)] or [19, Fact 3.1]) shows that 
/* — (•, z) is coercive. Consequently, all sublevel sets of tp are bounded. Hence, the 
sequence {\7 f {x'')} ken is bounded. By hypothesis, V/* is bounded on bounded 
subsets of V/(domA) and, according to (|1.7p . {a^'^ji.pRi is contained in V/(domA). 
Hence, the sequence x'^ — V/*(V/(x'^)), fc e N, is bounded. This proves Claim 1. 
Now we are going to prove the following: 

Claim 2: The sequence {^^^^en ^°^^ weak accumulation points and any such 
point is a solution of Problem 1.1. 

The space X being reflexive, there exists a weakly convergent subsequence {a;*'= }k£n 
of {x'^lfcgN- Let X =w-limfc^oo x'^'' ■ In order to show that S e C, denote y^ = a;*'= 
for each fc e N. According to (|4.22p . we have that y^ G Ci^, for every fc e N. 
By hypothesis w-lim/j^ooC'fe = C and this implies that x G C (see Definition 4.2). 
It remains to prove that 0* G Ax. To this end, observe that, according to (|2.10p . 
(|4?25| and ([^M]) we have 

(4.27) 0<D-f{z,u'')-Df{z,x''+'^) <Df{z,x'')-Df{z,x''+^), Vfc e N. 

As noted above, the sequence {Df{z, a;'°)}fegN is nonincreasing and nonnegative and, 
therefore, it converges. By (|4.27p . this implies that the sequence {D{z,u'^)}k^m 
converges and has the same limit as {Df{z,x'^)}keN- By (|4.23ll we also have that 

(4.28) Df{u^,x^)<Df{z,x^)~Df{z,u^), Vfc G N, 
and, thus, 

(4.29) \\m Df{u^,x^) = Q. 
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Since / is uniformly convex on bounded subsets of int doni /, it results that it is 
sequentially consistent too (cf. [36l Theorem 2.10]). Therefore, the equality ()4.29p 
implies that 

(4.30) lim ||M''-a;'=|| =0. 

Now, we distinguish two possible situations. First, suppose that condition (a) is 
satisfied. Note that, according to ()4.20p . we have 

C''' = V/(x'=)- V/(u'=), VfceN. 

Since V/ is uniformly continuous on bounded subsets of its domain, we deduce 
by (|4.30p that lim/j^oo ( — 0*. The operator A being sequentially weakly-strongly 
closed, this and the fact that {2:"'}^, j^ converges weakly to x imply that 0* £ Ax. 
Hence, a; is a solution of Problem 1.1 when condition (a) is satisfied. 

Alternatively, suppose that condition (6) is satisfied. Recall that, in this case, 
we necessarily have dom/ — X (cf. Remark 4.2(d)). By hypothesis (&), we have 
that C C dom A. By Assumption 1.1, we have that 

(V/ - A){C) C (V/ - A)(Co) C int dom/*, 

and, as shown above, x £ C. Hence, 

7^ (V/ - A){x) C int dom/* 

which implies that x G domA-f. From (|4.2p written with x^'' instead of y and u'* 
instead of u, we obtain that there exists u G A^x, such that 

(4.31) Dfiu'^u) KDfix'^x), VfceN. 
This implies 

(4.32) limmWfiu'^u) < liminfD/(a;'^S). 

k — ^00 k — 'oc 

Proposition 4.1, (|4.29p and (|4.32p imply that a; is a fixed point of A^ , that is, 
0* G Ax. Hence, x is a solution of Problem 1.1. This completes the proof of {{). 

The fact that if the Problem 1.1 has unique solution then the sequence {a;'^} , „ 
converges weakly to that solution is an immediate consequence of (i). Assume 
that the function / has sequentially weakly- weakly continuous gradient. We show 
next that, in this case, the sequence {a^'^/i.^fa can not have more than one weak 
accumulation point. Suppose by contradiction that this is not the case and that x' 
and x" are two different weak accumulation points of {a^'^}j,ppj ■ Let {a:*''}^pp^ and 
{■^"''° IfccM ^^ subsequences of {a;*'}, pjj converging weakly to x' and x" , respectively. 
By (i) combined with Lemma 3.3(&) it results that {x' , x"} C Sf{A, C) =:Nexp^(A). 
Hence, the inequality (|4.24p still holds for any z £ {x',x"}. It implies that the 
sequences {Df{x\x'')} ^ and {Df{x" ,x'')] ^ are convergent. Let a and b be 
their respective limits. For any fc S N we have 

Df{x',x^) ~ Df{x'\x^) = Df{x',x") + {Vfix") - Vfix''),x' - x") 

because of p.lOp . Replacing in this equation x^ by x^'' and letting k -^ 00 we 
deduce that 

(4.33) a-6 = L»/(x',x"), 
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because V/ is sequentially weakly-weakly continuous. A similar reasoning with x' 
and x" interchanged shows that 

h-a^Df(x",x'). 

Adding this equality with (|4.33|) we obtain that Df{x\ x") = 0. This can not happen 
unless x' = x" because the function / is strictly convex on C n int dom / as being 
Legendre. Thus, we reached a contradiction and this completes the proof of [ii). It 
is clear that [iii) follows from (i) and {ii) since the gradient of any convex function 
in a finite dimensional space is continuous on the interior of its domain (see, for 
instance, [STJ Proposition 2.8]). This completes the proof of the theorem. D 

4.3. Consequences of Theorem 4.1. If AT is a Hilbert space provided with the 
function / = ^||-|| , then Theorem 4.1 has a somewhat simpler form and even 
strong convergence of the sequence generated by the proximal-projection method 
can be sometimes ensured. Note that in this case the operator A : X ^ 2-^ is 
Df-coercive if and only if it is firmly nonexpansive in the sense that 

(4.34) eeArC^A, }^^^^^'--^)^ll^-^ll'- 

Clearly, if A has this property, then the operator A^ (which is exactly I — A) is 

f 
nonexpansive and, thus, Z?/-nonexpansive. Since in this situation Proj^ is exactly 

the metric projection operator Projcj., we obtain the following result: 

Corollary 4.1. Let X be a Hilbert space. Suppose that A : X -^ 2^ is a firmly 

nonexpansive operator (i.e., it satisfies (|4.34p ]. If C is a nonempty, closed and 

convex subset of X , if {Ck\)^fzm ''■^ '^ sequence of subsets of X satisfying (|4.18|) 

and such that Ck H dom^ is convex and closed and contains C for each A: e N, 

and if Problem 1.1 has at least one solution, then the sequence \X^^ u^ks generated 

according to the rule 

(4.35) a;° e Co n domA and x^^'^ e ^'^'^Ick+indomAi^^ - A^^), 

converges weakly to a solution of Problem 1.1. If int5/(A, C) 7^ 0, then ja;'^}, „ 
converges strongly. 

Proof. As noted above, the operator A satisfying (|4.34p is D/-coercive and Df- 
nonexpansive. Applying Theorem 4.1(6) with / = ^ ||-|| which has V/ = / (and, 
hence, has V/ sequentially weakly-weakly continuous) and taking into account that 
A^ = I — A, one deduces that the sequence {a^'^} 1,^™ converges weakly to a vector in 

Sf{A,C). The set Sf{A,C) = Nexpj-. A^ is convex and closed (cf. Lemma 3.1(a)). 
In the current circumstances, the inequality (|4.26p still holds for all z G Sf{A,C). 
It is equivalent to the condition 

||z-a;'=+i|| < ||z-x'=||, yz€Sf{A,C), Vfc G N. 

Therefore, one can apply Theorem 4.5.10 in [55] and this result implies that the 
sequence {a^'^jj.ppj converges strongly when int Sf {A, C) ^ 0. D 

It was pointed out in Subsection 3.3 that there is a strong connection between 
the proximal-projection method (|1.7p and the proximal point method (|3.23p - see 
Lemma 3.5. As far as we know, convergence of the proximal point method in reflex- 
ive Banach spaces was established for maximal monotone operators only. We use 
the connection between the proximal point method and the proximal-projection 
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method in order to obtain convergence of the proximal point method for opera- 
tors which are monotone with sequentially weakly-weakly closed graphs (but are 
not necessarily maximal monotone). Clearly, in spaces with finite dimension any 
monotone operator with closed graph and, in general, monotone operators with 
closed convex graphs have this property. The other requirement of the next corol- 
lary that dom_B/ should be convex is necessarily satisfied if V/* and V/ + B are 
surjective because, in this case, dom Bf = ran V/*o (V/ + i3) = X. This condition 
is sufficient without being necessary as the example preceding Lemma 3.5 shows. It 
can be easily verified that this also happens whenever Graph i? is convex and V/ 
is linear. Since the corollary is based on Theorem 4.1(a), the remarks preceding 
Theorem 4.2 concerning the implications of the hypothesis on the domains of / and 
/* still apply here. 

Corollary 4.2. Suppose that the following conditions are satisfied: 

(a) / is uniformly convex on bounded subsets of int dom /; 

(6) V/ is uniformly continuous on bounded subsets of int dom / as well as se- 
quentially weakly to weak continuous; 

(c) V/* is bounded on bounded subsets of int dom/*. 

If B : X ^r 2^ is a monotone operator with sequentially weakly-closed graph in 
X X X* , satisfying (|3.24p and such that dom Bf is convex, if B has at least one 
zero in int dom/, and if either 

(d) dom Bf is closed in X, 
or 

(e) V/ is bounded on bounded subsets of int dom /, 

then the sequences generated by the proximal point method ()3.23|) converge weakly 
to zeros of the operator B. 

Proof. We start by observing that, due to the boundedness on bounded subsets 
of its domain of V/* , we have that dom /* = X* - see Remark 4.2(d) . We first prove 
that the conclusion holds when dom. Bf is closed in X. Subsequently we will show 
that, if V/ is bounded on bounded subsets of int dom/, then dom Bf is necessarily 
closed in X and, thus, the conclusion is true in this case too. 

So, assume that dom Bf is closed in X. By Lemma 3.5, the proximal point 
method is identical to the proximal-projection method applied to the operator 
A[i?/] given by (|3.20p . Therefore, for proving the corollary in this case, it is suf- 
ficient to show that ^[B/] satisfies the requirements of Theorem 4.1(a). In order 
to do that it is sufficient to ensure that the operator A[i?/] is Z?/-coercive on its 
domain and has sequentially weakly-strongly closed graph. Observe that, accord- 
ing to Lemma 3.4, the operator A[Bf] is Z^z-coercive on its domain. It remains to 
show that A[Bf] is sequentially weakly-strongly closed. Let {y'^j^.^pj be a weakly 
convergent sequence contained in domA[_By] and denote ^ = A[Bf]y'^. Suppose 

that <(,> converges strongly in X* to some vector ^. Let y =w-limfe^oo y''- By 
L J fceN 

hypothesis (6), the sequence {^/(y'^)}^ j^ converges weakly in X* to ^ f{y)- Thus, 
the sequence 

(4.36) V/ {y^) -e= Wf o (V/ + B)-' o V/l {y'') 
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converges weakly in X* to V/(y) - ^. Denote u'' := V/* (v/ (y'') -<^'') and 
observe that, by (|4.36p . we have that 

(4.37) u' = [(V/ + B)-' o V/] (/■) = B//, Vfc e N. 

According to hypothesis (c), the sequence {u'''}j,gpj is bounded because the se- 
quence \ V/ (y*^) — f r" is bounded (as shown above this sequence is weakly 
I ^ ' J ken 

convergent). Let {w*'=}^ j^ be a weakly convergent subsequence of {w'^}^ p^j and 
let u be the weak limit of this subsequence. By (|4.37p we deduce that Vf{y^'') G 
(V/ + B) u^'' for all k E N, and thus we obtain 

(4.38) yfiy"") - yfiu'") e Bu'^, yk e N. 

By hypothesis (6), we have that 

(4.39) w- lim [Wfif") - Wf{u'-)] = Wf{y) - Wf{u). 

Since Graph i? is sequentially weakly- weakly closed and {«*'= }^pj^ converges weakly 
to u, the relations ^jT^ and ((i:^ imply that Vf{y) - V/(u) G Bu, i.e., Vf{y) e 
V/(u) + Su. Consequently, we have that 

(4.40) w=(V/ + B)-i(V/(y))=%y, 

because the operator Bf is single valued (cf. Lemma 3.4). On the other hand, by 
g36l), (I437|) . ([439]) and (14401) . we have that 

C = V/(2;) - w- hm Wfiu'") = V/(2/) - Wf(u) = V/(2;) - (V/ o Bf) y, 

h — *oo 

showing that {y,£,) G Graph yl[i3j^]. Hence, ^[i?/] is sequentially weakly-strongly 
closed and the proof, in this case, is complete. 

Now, assume that V/ is bounded on bounded subsets of int dom /. Then, by 
Remark 4.2((i), dom/ — X. We are going to show that, in this case, the set dom By- 
is closed. As shown above, if dom Bf is closed, then the conclusion holds. In order 
to prove that dom_B/ is closed, let l-^'^}, ^ be a sequence contained in dom Bf and 
converging in X to some vector z. Denote w'^ — Bfz^ . We claim that the sequence 
{w*^} is bounded. To show that, note that, since A[Bf] is _D^-coercive (cf. 

Lemma 3.4(c)), all solutions of Problem 1.1 are in the set Nexp-^^^^^,^ , {A[Bf]) 

(cf. Lemma 3.3). Taking into account that, by Lemma 3.4(a, 5), Bf = {A[Bf]Y and 
doTuBf = dom A[Bf\, we deduce that all solutions of Problem 1.1 are contained in 
Nexp^^j^^ Bf. Let z be such a solution. Then 

Df{z, w'') + Dfiw'', z^) < Df{z, z''), Vfc e N. 

According to the definition of the modulus of total convexity of / on the bounded 
set {z'^lj, , denoted '^/{{z'^} ^.^jq r)j ~ see [35] - we obtain 

(4.41) < i^f {{z'}^^^^ ; \\w' - zi) < Dfiw\z'^) < Df{z, z'^), Vfc e N, 

Since V/ is bounded on bounded subsets of X, the function / is also bounded 
on bounded subsets of X. Consequently, taking into account (|2.10p . we deduce 
that the sequence {Z3/(z,z'^)}, p, is bounded. Let M be an upper bound of this 
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sequence. Suppose by contradiction that the sequence {'u^''}i.piM contains a sub- 
sequence {w-'''}, p, such that hmfc_»oo ||w"'''|| = oo. Then there exists a positive 
integer fco such that for all integers fc > fco we have ||w-"= — z^'=|| > 1. By [551 
Proposition 2.1(m)] and (|4.4ip . we deduce that for any k > ko we have 

(4.42) 

The function / is, by hypothesis, uniformly convex on bounded subsets of int doni / 
and, consequently, it is also totally convex on bounded subsets of int dom / - cf. [361 

Theorem 2.10]. Therefore, vf I {^'^Itpia ; l) > 0. Taking this into account together 
with the fact that {z-"°}i.pRj is bounded (as being convergent) and letting k ^ oo 
in (|4.42p . we reach a contradiction. Hence, the sequence {w'^jj.gjj is bounded. Let 
{^''''ifceN ^^ ^ weakly convergent subsequence of {''w'^lijpRi and let w be the weak 
limit of this subsequence. According to (|3.19p . we have that 

(4.43) V/(z'^') - V/(w'=) e Bw'', Vfc G N. 

Since V/ is sequentially weakly-weakly continuous, and since B has sequentially 
weakly-weakly closed graph, the relation (|4.43p . written with Sk instead of /c, implies 
that V/(z) — V/(iD) e Bw. This shows that w — BfZ, that is, z E domBf. Hence, 
dom By is closed and the proof of the corollary is complete. D 

Another result which follows from Theorem 4.1 concerns a method of regularizing 
and solving classical variational inequalities in the form (|2.18p . Since the problem 
of solving ()2.18p may be ill-posed (in the sense that it may not have solutions or 
it may have multiple solutions) and, therefore, many algorithms for approximating 
solutions may not converge, or may converge only subsequentially, to solutions of 
the problem, one " regularizes" the original problem by solving an auxiliary problem 
which has unique solution and whose solution is in the vicinity of the solution set 
of the original problem, provided that the later is not empty. A regularization 
technique, which originates in the works of Tikhonov [72 and Browder [53], [27] . 
consists of replacing the original variational inequality (|2.18p by the regularized 
variational inequality 

(4.44) Find x E C n int dom / such that 

3^eBx: [{i + a\7f{x),y-x) > 0, Vy G Cndom/] , 

for some real number a > 0. If B is a monotone operator, then B + aV f is 
strictly monotone and, therefore, the variational inequality (|4.44p can not have 
more then one solution. Moreover, in many practically interesting situations, the 
variational inequality (|4.44p has solution even if the original variational inequality 
(I2.18P has not and, if a is sufficiently small, then the solution of (|4.44p is close 
to the solution set of the unperturbed variational inequality (j2.18p whenever the 
later has solutions. This is, for instance, the case (cf. [3 Theorem 3.2]) when the 
Banach space X is simultaneously uniformly convex and uniformly smooth and 
endowed with the Legendre function / :— - \\-\\^ for some p > 1 and B is maximal 
monotone. Theorem 4.1 allows us to prove the next corollary which extends the 
applicability of this regularization technique to reflexive Banach spaces which are 
not necessarily uniformly convex and uniformly smooth and to produce a weakly 
convergent algorithm for solving (|4.44p in this, more general setting, even if B is 
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not maximal monotone. This is of interest because closedness of the solution of 
()4.44p with small a > to the (presumed nonempty) solution set of the original 
variational inequality (|2.18p can be guaranteed even if B is not maximal monotone 
(cf. J, Theorem 2.1]). 

Corollary 4.3. Let B : X ^ 2^ be a monotone operator and let C be a 
nonempty, convex and closed subset of domi? fl int dom/. Suppose that f is uni- 
formly convex on bounded subsets of intdom/, V/* is bounded on bounded subsets 
of intdom/* and that, for some real number a > 0, we have that 

(4.45) 0^ ((l-a)V/-B)(C) C intdom/*, 

and the operator Proj^ [(1 — a)V/ — B)] is Df-firm. If one of the following con- 
ditions is satisfied: 

(a) X has finite dimension, dom f = X, V/ is uniformly continuous on bounded 
subsets of X, and B has closed graph and is bounded on bounded subsets of its 
domain; 

(b) V/ is bounded on bounded subsets of intdom/, the operator 
Proj^ [(1 — a)V/ — B)] is D j -nonexpansive; 

then the iterative procedure defined by 

(4.46) a;° e C and x''+^ e Proj^ [(1 - a)V/(x'=) - Bx''] , Vfc G N, 

is well defined and converges weakly to the necessarily unique solution of the vari- 
ational inequality ()4.44p . provided that such a solution exists. 

Proof. Well definedness of the procedure results from (|4.45p . The variational 
inequality (|4.44p can not have more than one solution since the operator B' := 
B + aV f is strictly monotone. According to Lemma 2.4, finding a solution of 
()4.44p is equivalent to finding a zero for the operator V :— V[B'\ C; f] defined by 
(f2T9l) . Note that 

(4.47) F = V/ - V/ o Proj^((l - a)V/ - B), 
and 

(4.48) T/-^=Proj^((l-a)V/~B) 
and that, by (|2.2ip applied to B' instead of B, we have 

(4.49) Proj^ (V/ -V)= Proj^((l - a)V f - B). 
Therefore, we can equivalently re-write the procedure (|4.46p as 

(4.50) x°eC and x''+^ e Pioi{j{Wf {x'') - Vx''), Vfc G N. 



This is exactly (|1.7p applied to V instead of A with the sequence of sets Ck = C 
for all k £ N. Now, suppose that condition (a) of our corollary is satisfied. In 
this case, if we show that the graph of V is closed m X x X and that V is Df- 
coercive, then Theorem 4.1(a) applies and leads to the conclusion of the corollary. 
Also, Theorem 4.1(6) implies that, if condition (&) of the corollary holds, then the 
procedure (|4.50p is weakly convergent to the unique solution of (|4.44p . provided 
that V is _D/-coercive. L)/-coercivity of V results in both cases from Lemma 3.2 
combined with (|4.48p and with our hypothesis that V^ =Proj^ [(1 — a)y f — B)\ 
is I?/-firm. So, it remains to prove that, under assumption (a) of the corollary, the 
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graph of V is closed. To this end, let {y'^j^.gfj be a sequence in domy and assume 

that this sequence converges to y e X. Let IS,'']- be the sequence 

L" ) ken 

e' = V/(/) - V/ o Proj^((l - a)V/(?/) - C'), 

where C & By'' for all fc G N. Suppose that linife^oo C = C- Then, by Lemma 2.1, 
we have 

v/(/-) - e' - fv/ o (v/ + Ncr'] ((1 - a)v/(/) - c'). 



VfceN. 



This implies that 

(4.51) (1 - a)V/(/) - C' e (V/ + Nc) [v/* {wf{y'') - e) 

Since B is bounded on the bounded set {y'^}i,pRii it follows that the sequence 

IC (■ is bounded. Let {C ] , ^m be a convergent subsequence of -^ C r and 
L J fcGN ... . . . '- J fe.eN 

let (^ be its limit. Since V/ and V/* are continuous on their respective domains and 

the normality operator Nc is maximal monotone (and, hence, has closed graph), 

(|4.5ip implies 

(1 - a)\/f{y) - C e (V/ + Nc) [V/* (V/(y) - $)] . 
Therefore, we have 

V/* (V/(y) -0= Proj^ [(1 - a)V/(y) - C] 
showing that 

C = V/(2/) - V/ o Proj^ [(1 - a)V/(2/) - C] G Vy. 

This completes the proof. D 

The requirement made in Corollary 4.3 that the operator V^ = 
Projp [(1 — a)V/ — B] should be D/-firm may seem unusual. Here are several 
examples which show that this condition is quite often satisfied without excessively 
costly demands on the operator B or the function /. The next example shows that 
Corollary 4.3(6) is applicable to solve variational inequalities in the form (I4.44[) 
when X is any Hilbert space, / = ^ ||-|| , C ^ domi? ^ X and S is a monotone 
operator which is either contractive with some constant 7 > or strongly monotone 
with some constant 5 > and even in more general conditions (when ()4.52p holds 
for some a € (0, ^) and (3 > 0). Obviously, in this setting, (|2.18p is exactly the 
problem of finding a zero of B. To follow the considerations in the examples below 
one should first note that by replacing in (|2.18p the operator B by (3B, where /3 is 
a positive constant, one obtains a variational inequality which is equivalent to the 
original one. 

Example 4.2. Let X be a Hilbert space and let f — ^ \\-\\ . Suppose that C = 
dom_B ~ X. If B is contractive with some constant 7 > or strongly monotone 
with some constant S > 0, then the operator V^ = Proj^ [(1 — a)V/ — /3B] is Df- 
firm and D f -nonexpansive whenever a G (0; 5) "^^^^^ (^ — \/oi{l^^^a}j~^ in the 
contractive case and /3 = (1 — 2a)S in the strongly monotone case. 

In order to prove this observe that, in the current setting, V/ = Proj^ = /. 
Also, the Z?/-firmness condition (|3.14p for T = V^ is equivalent to (|4.34p and this 
is exactly 

([(1 - a)x - P^] - [(1 - a)y -Pr,],x-y) 
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l2 



>\\[[i^a)x-m-[{^-o^)y~H\\' 

which can be cquivalently re-written as 

(4.52) a{l~a)\\x-v\\^ + {l-2a)P{£,-T^,x-y)>0'\\i-T^\\ 

for all pairs (x,^), (y,?]) G Graphi?. Moreover, if V^ is Z?/-firni, then it is also Dj- 
nonexpansive. Due to the monotonicity of B the inequality (|4.52p holds whenever 
a € (0, 5) and 

(4.53) a{l~a)\\x~y\\>l3\\S,-ri\\, V(a;,C) e Graph S,V(y, 77) G Graph B. 
If B is contractive with constant 7, then we have 

l\\x- v\\ > lie - 7?ll , V(x,e) e Graph S,V(y, 77) e GraphS 



and multiplying this inequality by /3 = \J a{\ — 0)7 ^ we deduce that (|4.53p holds. 
Similarly, note that (|4.52p is satisfied when 

(4.54) {l-2a){i-f^,x~y) >/3||C-?7l|^ V(a;,0 e Graph B, V(2/, ry) e GraphB. 
If B is strongly monotone with constant (5, then 

(e-77,a;-y) > (5||C-?7f , V(a;,^) G GraphS, V(y, 77) G GraphB. 
Multiplying this inequality by (1 — 2a) we deduce that (|4.54p holds in this case. D 

In situations in which X is not a Hilbert space, or X is a Hilbert space but the 
monotone operator B does not satisfy (|4.52|) for some a G (0, ^) and /3 > 0, the 
question of how to choose the "regularization function" / and the " regularization 
parameter" a in the perturbed variational inequality (I4.44[) in order to force D f- 
firmness and/or _D^-nonexpansivity on V^ = Projp((l — a)V/ — B) is relevant. 
Here are examples of situations when V^ is L'j-firm even if X is not a Hilbert 
space. 

Example 4.3. // the monotone operator B : X ^ 2^ and the nonempty 
closed convex subset C of domi? n intdom/ have the property ()4.45p . then V^ = 
Projp((l — a)V/ — B) is Df-firm in any of the following situations: 

(a) C = X, a G (0, 1) and the operator aV f + B is Df -coercive on its domain; 

(b) X is a Hilbert space, C = X, / = 2 ll'll ; o; G (0; 5] ^^'^ ^ ^^ contractive with 
constant a(l — a). 

In order to show this observe that in case (a), according to Lemma 2.1, we have 
V^ — {aWf + B)f and, by Lemma 3.2(c), the conclusion follows. In case (&) an 
easy calculation shows that the operator aV/ + B is _D/-coercive and, therefore, 
the conclusion of case (a) applies. D 

Example 4.3(a), in conjunction with Corollary 4.3, leads to an algorithm for 
solving the variational inequality (j4.44p whenever is possible to find a Legendre 
function / and a number a G (0, 1) such that aV/ + _B is _D/-coercive on its 
domain. Example 4.3(6) covers the class of expansive operators in Hilbert spaces 
for which, as far as we know, few methods of finding zeros are available. 

Example 4.4. If the monotone operator B : X ^ 2^ and the nonempty closed 
convex subset C of dom B D int dom / have the property (I4.45P and if 

(4.55) ([aV/(x)+e]-[aV/(y)+77], 

Proj^ [(1 - a)V/(x) - e] - Proj^ [(1 " a)Vfiy) - v]) > 0, 
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for any (x, £,) and (y, rj) in Graph _B, then the operator V^ ~ Projp [(1 — a)V/ — B] 
is Df-firm. In particular, this happens in any of the following situations: 

(a) B=[l~a)Vf and a e (0, i); 

(b) The operator P : X x X* -^ 2-^'^^, defined by 

(4.56) F(z, C) = (0*, Proj^((l - a)V/(z) - 0), 

for some a G (0, 1) is monotone when X x X* is provided with the norm ||(z, C,)\\ = 

1 /? 
(Ikll +IICIL) o'^'^ with the duality pairing /(z, C), (C', z')^ = {C 1 ^} + {Cj^') 

(and, therefore, its dual is isometric with X* x X). 
Let X, 1/ e dom V^ and let ^ e Bx and rj e Sy. Denote 

x' = Proj^((l - a)Vf{x) ~ S) and y' = Proj^((l - a)V/(zj) - 77). 

By Lemma 2.1 we have 

(1 - a)Vf{x) - e - Vf{x') e Nc{x') and (1 - a)V/(y) - r; - V/(y') G 7Vc(y'), 

which imply 

((1 - a)V/(x) - ^, y' - x') < {Vf{x'),y' - x') 
and, respectively, 

((1 - a)V/(y) - r,, a;' - y') < (V/(y'), a^' " v') ■ 
Summing up the last two inequalities we obtain that 

(4.57) ([(1 - a)Vf{x) ~i]~ [(1 - a)V/(y) ~ii\,x'- y') 

>(V/(x')-V/(y'),x'-y'). 
The Z?/-firmness condition (|3.14p for the operator T = V^ is exactly 

(V/(x) - V/(y), x' - y') > (V/(a;') - V/(y'), x' - y') . 

According to (|4.57|) . this is satisfied when 

{^f{x) - V/(y), x' - y') > ([(1 - a)Vf{x) - C] - [(1 - «)V/(y) - r;] , x' - y'> 

and this last inequality is equivalent to (|4.55[) . Hence, if (|4.55p holds, then the 
operator V^ is ZJ^-firm. In case (a) we have that aSJ f + B = (1 — a)Wf — B and 
using the monotonicity of Proj^ (cf. [311 Theorem 4.6]) one deduces that (|4.55p 
holds. Suppose that we are in case (d) and the operator P, given by (|4.56p . is 
monotone. Then observe that 

([aV/(x)+C]-[aV/(y) + 77], 
Proj^ [(1 - a)Vfix) - C] - Proj^ [(1 - a)V/(y) - r,]) = 

{{x, aVf{x) +0- (y, «V/(y) + 7?), P(x, - P(y, v)) , 
and that the last expression is nonncgativc due to the monotonicity of P (see |36l 
Proposition 4.7]). Hence, (|4.55p holds in this case too. D 

Note that problem (|4.44p in which B — j3V f for some /3 > is equivalent to 
the problem of finding the minimizer of / over C. The facts observed in Example 
4.4(a), in conjunction with Corollary 4.3, leads to a proximal-projection method of 
finding that minimizer, provided that / satisfies the other requirements there. Ob- 
viously, the effectiveness of that method, as well as of the other methods discussed 
in this work, depends on the possibility of computing proximal projections onto 
C. Algorithms for computing proximal projections are presented in [5], [H] and 
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[36] . Compared with already classical projection methods (see p?) in which the 
iterations are, usually, of the form x'^^^ = Projp(a;'^ — AfcV/(x'^)) with a converging 
to zero positive step size Afc, the proximal-projection method presents the advan- 
tage of not requiring arbitrarily small step sizes which, in practical applications, 
may force the procedure to became stationary long before the iterates are close to 
the minimizer of / (due to the computer identification of AfeV/(a;'^) with the null 
vector). 

In general, it would be nice to have a Legendre function / for which the condition 
in Example 4.4(6) is satisfied. Whether such a function exists in a Banach space 
X is an open question. This question is relevant because such a function /, if any, 
would be a "universal regularizer" for variational inequalities in the form ()2.18p in 
the sense that it would be such that the regularized variational inequality ()4.44p 
will be solvable by the proximal-projection method, no matter how the monotone 
operator B is, provided that (I4.44p has solution. 

5. Convergence of the Proximal-Projection Method in Spaces of 

Finite Dimension 

Theorem 4.1 and its corollaries ensure weak and, sometimes, strong convergence 
of the proximal-projection method to solutions of the Problem 1.1 under conditions 
which, besides the Z3/-coercivity of the operator A, require sequential weak-strong 
closedness of the Graph A or, alternatively, D/-nonexpansivity oi A-^ . In this section 
we show that, when the space X has finite dimension, some of these requirements 
can be dropped or weakened. This is possible due to the validity in spaces of finite 
dimension of another generalization of Opial's Lemma which we present below. 

5.1. Another variant of a generalized Opial's Lemma. The following result 
applies to operators T : X ~* 2^ which are not necessarily ZJ^-nonexpansive, but 
satisfy condition (|5.2p below which is more general than Z3/-firmness (compare con- 
dition (|5.2p with Definition 3.3). It is interesting to observe that, if / is uniformly 
convex on bounded subsets of int dom / and if T has closed graph, then the conclu- 
sion of the next result holds even if the hypothesis that u satisfies (|5.2p is removed. 
This happens because the equality in (|5.ip implies that the sequences \Z } , ,, 
and {m'^}^ j^ converge to the same limit z and, then, closedness of the graph of T 
guarantees that z £Tz. 

Proposition 5.1. Suppose that the space X has finite dimension, f is uniformly 
convex on bounded subsets of int dom / and T : X —> 2^ is an operator satisfying 
condition (I3.13p . Let {z^}k£fi be a sequence in domT converging to an element 
z G domT. If for some sequence {w'^jfegN satisfying 

(5.1) (Vfc e N : u^ e Tz^) and lim Df{u^,z^) = 0, 



k- 



there exists u € Tz such that 



(5.2) liminf {Vf{u'') - Wf{u), u'' - u) < 0, 

k — >oo 

then the vector z is a fixed point of T . 

Proof. Since the function / is convex and differentiable on int dom /, the gra- 
dient V/ is continuous on int dom/. This fact and the strict convexity of / on 
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int dom / imply that 

(5.3) lim (V/(z'=) - Vfix), z'= - x) = (V/(z) - V/(a;), z - x) > 0. 

k — *oo 

whenever x G (int dom /) \ {z} . The function / being uniformly convex on bounded 
subsets of int dom/, it is also sequentially consistent (cf. [3S1 Theorem 2.10]). 
Therefore, the equality in (|5.ip implies that 

(5.4) lim |U''-u''|| =0. 

k — ^oo 

Consequently, the sequences {^'^jfegN and {w'^j^gN converge to the same limit z. By 
condition p.l3p . the boundcdncss of {u'^j/cgpf and the continuity of V/ we deduce 
that 

(5.5) lim (V/(z'=)-V/(z),z'=-z) =0= lim (V/(z'^) - V/(z),u'= - u). 

Note that 

(5.6) {Vf{u'^)-Vfiu),u''-u) = {Vf{u'')-\/fiu),u'^-z'^) 

+ {Vf{u'')-S/f{u),z'^-u) 

= {\/f{u')-\/f{u),u''-z'') 

+ {Wf{z'')-Wf{u),z'^-u) 

+ (V/(m'=)-V/(z'=),z'=-«). 

Since the sequence {'^f{u'^)}keN is bounded, it follows from (|5.4p that the first 
term of the last sum in ()5.6p converges to zero. The second term of the same sum 
is nonnegative because of the monotonicity of V/. Taking the limit as fc — j oo on 
both sides of (|5.6p . we obtain that 

(5.7) lim (V/(u'=) - V/(u), u'' - u) > lim {Wf{z'') ~ V/(u), z'^^u). 

k — ^oo k — >oo 

In order to conclude the proof, suppose by contradiction that z ^ T{z). Then, 
u^ z and, therefore, by (|5.5p . (j5.2p . (|5.7p and (|5.3p . respectively, we obtain 

= lim (V/(z'=)-V/(z),z'^-z)= lim (V/(z''VV/(z),u'=-u) 

k^oQ k — !-oo 

> lim (V/(u'=) - Wf{u),u'' - u) > lim {Wf{z'') - V/(u), z'' ~ u) > 0, 

k — *C30 k — ^oo 

which is a contradiction. D 

5.2. A convergence theorem for the proximal-projection method in spaces 
of finite dimension. The following theorem shows that, in finite dimensional 
spaces, convergence of the proximal-projection method to solutions of Problem 1.1 
can be ensured with lesser requirements on the operator A in addition to the Dj- 
coercivity than those involved in Theorem 4.1 and its corollaries. 

Theorem 5.1. Suppose that the space X has finite dimension, f is uniformly 
convex on bounded subsets of int dom/, V/* is bounded on bounded subsets of 
V/(domA) and that ()1.6p . Assumption 1.1 hold and 

(5.8) C n dom A = w-Iimfe^oo {Ck n dom A) . 



If the Problem 1.1 has at least one solution, if the operator A : X — > 2-^ is D 



f- 



coercive on 

k=0 



Q — \ j Ck and if C O dom A is closed, then the sequences generated 
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by the proximal-projection method (|1.7|) are well defined and converge to solutions 
of the Problem 1.1. 

Proof. Well definedness of the sequences generated by (|1.7p follows from (|1.6p 
and Assumption 1.1. Suppose that, for each fc S N, C and u^ are as m ()4.19|) 
and (I4.20p . respectively. Then, clearly, condition ()4.21|) holds too. The operator 
A being Z)^-coercive on its domain, the operator A^ is Z?/-firm (cf. Lemma 3.2). 
This means that 

(5.9) {^fiu'') - Vf{u)y ~u)< (V/(x^-) - V/(x), «^ - u) , 

for any pair (x,u) G Graph A-' and for any fc G N. Now, repeating without change 
the arguments in the proof of Theorem 4.1 one can see that Claim 1 proven there 
still holds in our setting and implies that the sequence {a^'^j^.^j^ is bounded. Let 
{cc*'-}^ j^ be a convergent subsequence of {a^'^j^.gpj and let x be its limit. An 
argument identical to that made in the proof of Theorem 4.1 (Claim 2) for the 
same purpose shows that x G C and (|4.30p holds. According to (|5.8p . since a;*'= G 
Ci^ n domA, it results that a; G C n dom^. Hence, Ax ^ 0. Writing (|5.9p for 
ik instead of k and x instead of x, and for any u G Ax, letting in the resulting 
inequality k —t oo and taking into account that V/ is continuous on intdom/, we 
deduce that 

(5.10) liminf(V/(u"=)- V/(u),u"' -u) < 0. 

k — ^CxD 

This shows that the sequence {w"°}j,pj>[ satisfies (|5.2p ioi T — A^ , z = x and 
any u G A^ x. Also, by (|4.30p and (|5.10p . the sequences {2;*''}j.ppj and {u*''}^, 
satisiy (|5.ip . Hence, Proposition 5.1 applies to T = ^^, z — x and u G A^ x. By 
consequence, we have that S is a fixed point of A^ and, hence, a zero of A. It 
remains to prove that x is the only accumulation point of the sequence {2;'^}j.ppj . 
The proof in this respect reproduces without modifications the arguments made in 
the proof of Theorem 4.1 in order to show that the sequence {a;*^}^ pj has a single 
weak accumulation point when V/ is sequentially weakly-weakly continuous. D 

5.3. Consequences of Theorem 5.1. Using Theorem 5.1 instead of Theorem 4.1 
we can prove again Corollary 4.2 and Corollary 4.3 in a finite dimensional setting 
with different and less demanding conditions on A. Here is the new version of 
Corollary 4.2. 

Corollary 5.1. Suppose that the space X has finite dimension, f is uniformly 
convex on bounded subsets of intdom/ and V/* is bounded on bounded subsets of 
intdom/*. If B : X ^ 2^ is a monotone operator satisfying (|3.24p and having at 
least one zero, and if any of the following conditions holds 

(a) ran (V/ + B) is closed in X* and V/* (ran (V/ + B)) is convex; 

(b) ran {V f + B) = X* ; 

then the sequences generated by the proximal point method (|3.23p converge to zeros 
of the operator B. 

Proof. Recall that in this setting dom/* = X* (cf Remark 4.2(d)). Observe 
that, according to Lemma 3.4, we have that 

dova A[Bf] = dom 5/ = V/* (ran (V/ + B)) . 

Therefore, if condition (a) holds, the set dom A [Bf] — dovn Bf is closed and convex. 
Since condition (6) implies (a), this remains true when (6) holds. Again by Lemma 
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3.4, the operator A[Bf] is I?/-coercive on its domain. Consequently, the operator 
A :— A [Bf] satisfies the requirements of Theorem 5.1 when C ^ Ck ~ X for all 
k E N. Applying Theorem 5.1 to A[i?/] and taking into account Lemma 3.5 the 
conclusion follows. D 

Now we give another variant of Corollary 4.3(a) in which the condition that A 
should have closed graph is replaced by less demanding requirements. 

Corollary 5.2. Let B : X ^ 2^ be a monotone operator. Suppose that 
the space X has finite dimension, f is uniformly convex on bounded subsets of 
intdom/ and V/* is bounded on bounded subsets of intdom/*. If C is a closed 
convex subset of dom_B n intdom/ such that, for some real number a > 0, 

(5.11) 7^ ((1 - a)Vf - B){C) C intdom/*, 

and the operator Proj^ o [(1 — q;)V/ — B)] is Df-firm, then the iterative procedure 
given by (|4.46p is well defined and converges to the necessarily unique solution of 
the variational inequality (j4.44p . provided that such a solution exists. 

Proof. Since (|4.48p and (|4.49p still hold, the operator V given by (|4.47p is 
Z^Z-coercive on its domain (cf. Lemma 3.4). By (|5.1ip and by the fact that C C 
domB n intdom/, it results that C C domV. Hence, Theorem 5.1 applies to the 
operator A — V and the sets Ck = C and the conclusion follows. D 
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